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Abstract. In this paper, we study spectral properties of the one dimensional periodic Schrodinger 
operator with an adiabatic quasi-periodic perturbation. We show that in certain energy regions the 
perturbation leads to resonance effects related to the ones observed in the problem of two resonating 
quantum wells. These effects affect both the geometry and the nature of the spectrum. In particular, 
they can lead to the intertwining of sequences of intervals containing absolutely continuous spectrum 
and intervals containing singular spectrum. Moreover, in regions where all of the spectrum is expected 
to be singular, these effects typically give rise to exponentially small "islands" of absolutely continuous 
j spectrum. 
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Resume. Cet article est consacre a l'etude du spectre d'une famille d'operateurs quasi-periodiques 
obtenus comme perturbations adiabatiques d'un operateur periodique fixe. Nous montrons que, dans 
certaines regions d'energies, la perturbation entrame des phenomenes de resonance similaires a ceux 
observes dans le cas de deux puits. Ces effets s'observent autant sur la geometrie du spectre que sur sa 

■ nature. En particulier, on peut observer un entrelacement de type spectraux i.e. une alternance entre 
du spectre singulier et du spectre absolument continu. Un autre phenomene observe est l'apparition 
d'ilots de spectre absolument continu dans du spectre singulier dus aux resonances. 
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00 ■ 0. Introduction 

O 

t^- | The present paper is devoted to the analysis of the family of one-dimensional quasi-periodic 

■ Schrodinger operators acting on L 2 (R) defined by 

d 2 

Oh' (0.1) H z , e = -— — + V(x - z) + a cos (ex). 



dx 2 

We assume that 

(HI): V : R — ► R is a non constant, locally square integrable, 1-periodic function; 
(H2): e is a small positive number chosen such that 2ir/e be irrational; 
^ ' (H3): z is a real parameter; 

^ . (H4): a is a strictly positive parameter that we will keep fixed in most of the paper. 

As e is small, the operator (0.1) is a slow perturbation of the periodic Schrodinger operator 
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(0.2) H„ = ~ + V(x) 

acting on L 2 (R). To study (0.1), we use the asymptotic method for slow perturbations of one- 
dimensional periodic equations developed in [11] and [9]. 

The results of the present paper are follow-ups on those obtained in [12, 8, 10] for the family (0.1). 
In these papers, we have seen that the spectral properties of H Z)£ at energy E depend crucially on 
the position of the spectral window J-{E) := [E — a, E + a] with respect to the spectrum of the un- 
perturbed operator Hq. Note that the size of the window is equal to the amplitude of the adiabatic 
perturbation. In the present paper, the relative position is described in figure 1 i.e., we assume that 
there exists J, an interval of energies, such that, for all E G J, the spectral window J~(E) covers the 
edges of two neighboring spectral bands of Hq (see assumption (TIBM)). In this case, one can say 
that the spectrum in J is determined by the interaction of the neighboring spectral bands induced by 
the adiabatic perturbation. 
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The central object of our study is the monodromy equation, a finite difference equation determined 
by the monodromy matrix for the family (0.1) of almost periodic operators. The monodromy matrix 
for almost periodic equations with two frequencies was introduced in [12]. The passage from (0.1) to 
the monodromy equation is a non trivial generalization of the monodromization idea from the study 
of difference equations with periodic coefficients on the real line, see [3]. 
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Let us now briefly describe our results and the heuristics underlying them. Let E 
E(k) be the dispersion relation associated to Hq (see section 1.1.2) ; consider 
the real and complex iso-energy curves, respectively and T, defined by 

(0.3) r M := {(C, k) G M 2 ; E(k) + a ■ cos(C) = E}, 

(0.4) T := {(£ k) G C 2 ; E(«) + a ■ cos(C) = E}. 

The dispersion relation k h- > E(«) being multi- valued, in (0.4), we ask that the 
equation be satisfied at least for one of the possible values of E(«). 
The curves T and are both 27r-periodic in the k- and ("-directions; they are 
described in details in section 10.6. The connected components of T]r are called real branches of T. 

Consider an interval J such that, for E G J, the 
assumption on the relative position of the spectral win- 
dow and the spectrum of Hq described above is satis- 
fied (see figure 1). Then, the curve Tk consists of a 
infinite union of connected components, each of which 
is homeomorphic to a torus ; there are exactly two 
such components in each periodicity cell, see figure 2. 
In this figure, each square represents a periodicity cell. 
The connected components of Tk are represented by 
full lines; we denote two of them by 70 and 7^. 
The dashed lines represent loops on T that connect 
certain connected components of Tj&; one can distin- 
guish between the "horizontal" loops and the "ver- 
tical" loops. There are two special horizontal loops 
denoted by j hfi and j h y, the loop j h ,o (resp. 7^) 
connects 70 to 7^ — (2-7T, 0) (resp. 70 to In the 

same way, there are two special vertical loops denoted 

by 7^,0 an d 7i,,7r ; the loop 7^0 (resp. 7 Uj7r ) connects 70 to 70 + (0, 2n) (resp. 7^ to 7^ + (0, 2ix)). 

The standard semi-classical heuristic suggests the following spectral behavior. To each of the 
loops 70 and 7^, one associates a phase obtained by integrating the fundamental 1-form on T along 
the given loop; let <J>o = $o(E) (resp. & n = $ n (E)) be one half of the phase corresponding to 70 
(resp. 77,-). Each of these phases defines a quantization condition 




Figure 2: The adiabatic phase space 
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Each of these conditions defines a sequence of energies in J, say (Eq)i and (E% ^)i>- For e sufficiently 
small, the spectrum of E z s in J should then be located in a neighborhood of these energies. 

Moreover, to each of the complex loops 7/^0, 7/i l7r , 7^,0 and 7^,77, one naturally associates an action 
obtained by integrating the fundamental 1-form on V along the loop. For b G {0,7r} and a G {v,h}, 
we denote by S a b the action associated to 7^ multiplied by i/2. For J?£l, all these actions are real. 
One orients the loops so that they all be positive. Finally, we define tunneling coefficients as 

ta,b = e~ Sa ' b/e , be {0,vr}, ae{v,h}. 

When the real iso-energy curve consists in a single torus per periodicity cell (see [12]), the spectrum 
of H z<e is contained in a sequence of intervals described as follows: 

• each interval is neighboring a solution of the quantization condition; 

• the length of the interval is of order the largest tunneling coefficient associated to the loop; 

• the nature of the spectrum is determined by the ratio of the vertical tunneling coefficient to 
the horizontal one: 

— if this ratio is large, the spectrum is singular; 

— if the ratio is small, the spectrum is absolutely continuous. 
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In the present case, one must moreover take into account the possible interactions between the tori 
living in the same periodicity cell. Similarly to what happens in the standard "double well" case 
(see [14, 24, 15]), this effect only plays an important role when the two energies, generated each by 
one of the tori, are sufficiently close to each other. In this paper, we do not consider the case when 
these energies are "resonant", i.e. coincide or are "too close" to one another, but we can "go" up to 
the case of exponentially close energies. 

Let Eq be an energy satisfying the quantization condition (0.5) defined by §q; let S be the distance 
from Eq to the sequence of energies satisfying the quantization condition (0.5) defined by We now 
discuss the possible cases depending on this distance. Let us just add that, as the sequences of energies 
satisfying the quantization equation given by $o or ^tt play symmetric roles, in this discussion, the 
indexes and tt can be interchanged freely. 

First, we assume that, for some fixed n > 1, this distance is of order at least e n . In this case, 
near Eq, the states of the system don't "see" the other lattice of tori, those obtained by translation 
of the torus j n ; nor do they "feel" the associated tunneling coefficient t V)7T . Near Eq, everything is 
as if there was a single torus, namely a translate of 70, per periodicity cell. Near Eq, the spectrum 
of H Z E is located in a interval of length of order of the largest of the tunneling coefficients t v q and 
tfi = th,oth,-K (see section 1.3.3). And, the nature of the spectrum is determined by quotient t v $/th- 
So, in the energy region not too close to solutions to both quantization conditions in (0.5), we see 
that the spectrum is contained in two sequences of exponentially small intervals. For each sequence, 
the nature of the spectrum is obtained from comparing the vertical to the horizontal tunneling coef- 
ficient for the torus generating the sequence. As the tunneling coefficients for both tori are roughly 
"independent" (see section 1.7.5), it may happen that the spectrum for one of the interval sequences 
be singular while it be absolutely continuous for the other sequence. If this is the case, one obtains 
numerous Anderson transitions i.e., thresholds separating a.c. spectrum from singular spectrum (see 
figure 5(b)). 

Let us now assume that 5 is exponentially small, i.e. of order e~ ri ^ £ for some fixed positive r] (not 
too large, see section 1.6). This means that we approach the case of resonant energies. Note that, this 
implies that there is exactly one energy E n satisfying (0.5) for (Ev that is exponentially close to Eq; 
all other energies satisfying (0.5) for are at least at a distance of order e away from Eq. 
Then, one can observe two new phenomena. First, there is a repulsion of Iq and I n , the intervals 
corresponding to Eq and E n respectively containing spectrum. This phenomenon is similar to the 
splitting phenomenon observed in the double well problem (see [14, 24, 15]). Second, the interaction 
can change the nature of the spectrum: the spectrum that would be singular for intervals sufficiently 
distant from each other can become absolutely continuous when they are close to each other, see 
Fig. 5(a). To explain this phenomenon, assume, for simplicity, that t v ^ and t v>n , the "vertical" 
tunneling coefficients associated to the tori 70 and 7^, are of the same order (in e), i.e. t v ,o ~ t V{K ~ 
t v . Then, if \Eq — E n \ ~ e n , on each of the intervals Iq and /„-, the nature of the spectrum is 
determined by the same ratio t v /th- I£ \Eq — E n \ ~ e~ v ^ £ , the two arrays of tori begin to "feel" one 
another: they form an array for which the tori from both arrays play equivalent roles. In result, 
the "horizontal" tunneling becomes stronger: it appears that th has to be replaced by the effective 
"horizontal" tunneling coefficient th,eS — th / dist(EQ, E n ), and the ratio t v /th has to be replaced by 
tv/th,eS- So, the singular spectrum on the intervals Iq and I n "tends to turn" into absolutely continuous 
one. 

There is one more case that will not be discussed in the present paper: it is the case when 
5 ~ e^l^ with no restriction on rj positive or, even, when 5 vanishes. This is the case of strong 
resonances; it reveals interesting new spectral phenomena and is studied in detail in [7]. 

1. The results 

We now state our assumptions and results in a precise way. 

1.1. The periodic operator. This section is devoted to the description of elements of the spectral 
theory of one-dimensional periodic Schrodinger operator Hq that we need to present our results. For 
more details and proofs we refer to section 6 and to [6, 13]. 
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1.1.1. The spectrum of Hq. The spectrum of the operator Hq defined in (0.2) is a union of countably 
many intervals of the real axis, say [E 2n +i, E 2n +2\ for n £ N , such that 

El < E2 < E3 < £4 . . . E2n < E2n+l < E2 n +2 < • • • , 

E n — > +00, n — ► +00. 

This spectrum is purely absolutely continuous. The points (Ej)j^ are the eigenvalues of the self- 
adjoint operator obtained by considering the differential polynomial (0.2) acting in L 2 ([0,2]) with 
periodic boundary conditions (see [6]). The intervals [E 2n+ \, E 2n+2 ], n £ N, are the spectral bands, 
and the intervals (E 2n , -E^n+i); u6N*, the spectral gaps. When E 2n < E 2n +\, one says that the n-th 
gap is open; when [E 2n _i, E 2n ] is separated from the rest of the spectrum by open gaps, the n-th band 
is said to be isolated. 

From now on, to simplify the exposition, we suppose that 
(O): all the gaps of the spectrum of Hq are open. 

1.1.2. The Bloch quasi-momentum. Let x 1— > ip(x, E) be a non trivial solution to the periodic Schrodin- 
ger equation H^ip = Eijj such that, for some fx £ C, ^ (x + 1, E) = flip (x, E), \/x £ R. This solution 
is called a Bloch solution to the equation, and fi is the Floquet multiplier associated to tp. One may 
write u = exp(ifc); then, k is the Bloch quasi-momentum of the Bloch solution tj). 

It appears that the mapping E i— > k(E) is an analytic multi-valued function; its branch points are the 
points Ei, E2, E3, . . . , E n , .... They are all of "square root" type. 

The dispersion relation k •— ► E(/c) is the inverse of the Bloch quasi-momentum. We refer to section 6.1.2 
for more details on k. 

1.2. A "geometric" assumption on the energy region under study. Let us now describe the 
energy region where our study will be valid. 

The spectral window centered at E, T{E), is the range of the mapping £ £ R ^ E — acos(^). 

In the sequel, J always denotes a compact interval such that, for some n £ N* and for all E £ J, one 
has 

(TIBM): [E 2n ,E 2n+1 ] C F{E) and T{E) c]£ 2 n-i, E 2n+2 [. 
where J~{E) is the interior of J~{E) (see figure 1). 

Actually, in the analysis, one has to distinguish between the cases n odd and n even. From now on, 
we assume that, in the assumption (TIBM), n is even. The case n odd is dealt with in the same way. 
The spectral results are independent of whether n is even or odd. 

Remark 1.1. As all the spectral gaps of Hq are assumed to be open, as their length tends to at 
infinity, and, as the length of the spectral bands goes to infinity at infinity, it is clear that, for any non 
vanishing a, assumption (TIBM) is satisfied in any gap at a sufficiently high energy; it suffices that 
this gap be of length smaller than 2a. 

1.3. The definitions of the phase integrals and the tunneling coefficients. We now give 
precise definitions of the phase integrals and the tunneling coefficients appearing in the introduction. 

1.3.1. The complex momentum and its branch points. The phase integrals and the tunneling coeffi- 
cients are expressed in terms of integrals of the complex 
momentum. Fix E in J. The complex momentum C ^ n 
k(C) is defined by 

(1.1) k(C) = k(E- acos(C)). 

As k, k is analytic and multi- valued. The set T defined 
in (0.4) is the graph of the function k. As the branch 
points of k are the points (Ei)^^, the branch points of k 
satisfy 1 1 

(1.2) E- acos(C) = Ej, j £ N*. 
As E is real, the set of these points is symmetric with res- Figure 3: The branch points 

pect to the real axis, to the imaginary axis; it is 27r-periodic in £. All the branch points of k lie in the 
set arccos(R) which consists of the real axis and all the translates of the imaginary axis by a multiple 
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As the branch points of the Bloch quasi-momentum, the branch points of k are of "square root" type. 

Due to the symmetries, it suffices to describe the branch points in the half-strip {£; Im£ > 
0, < ReC < Ti"}- These branch points are described in detail in section 7.1.1. In figure 3, we show 
some of them. The points (Cj)j satisfy (1.2); one has 

< (2n < C2n+1 < TT, < Im ( 2 n+2 < Im C2n+3 < " " , < Im ( 2 n-l < • • • < Im Cl- 

1.3.2. The contours. To define the phases and the tunneling coefficients, we introduce some integration 
contours in the complex £-plane. 

These loops are shown in figure 3 and 4. The loops 70, %, 7/^0 > 7/i,7r> %,o and 7„ j7r are simple loops 
going once ar ound respec tively the intervals [~C2n,C2n], [C2n+i, 27r-C 2 n+i]. [-C271+1, -C,2n], [C2n, C2n+i], 
[C2n-1> C271-1] and [(2n+2 , C2n+2j ■ 

In section 10.1, we show that, on each of the above loops, one can fix a continuous branch of the 
complex momentum. 

Consider T, the complex iso-energy curve defined by (0.4). Define the projection II : (£, k) G T 1— ► £ G 
C. The fact that, on each of the loops 70, 7^, 7/i,o> lh,n, lv,o and t^tt, one can fix a continuous branch 
of the complex momentum implies that each of these loops is the projection on the complex plane of 
some loop in T i.e., for 7 G {%, j^, 7^0, 7/1,70 Tu,o 5 7«,7r}> there exists 7 C T such that 7 = 11(7). In 
sections 10.6.1 and 10.6.2, we give the precise definitions of the curves 70, 7^, 7h,o> 7/i,tt> 7u,o and 7„ j7r 
represented in figures 3 and 2 and show that they project onto the curves 70, j n , 7h,cb 7/i,7tj %,o and 
7u j7r respectively. 

1.3.3. The phase integrals, the action integrals and the tunneling coefficients. The results described 
below are proved in section 10. 



Let v € {0,7r}. To the loop 7^, we associate the phase 
integral Q u defined as 



(1.3) 
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where k is a branch of the complex momentum that is 
continuous on 7^. The function E \— > § U (E) is real ana- 
lytic and does not vanish on J. The loop 7^ is oriented 
so that $ U (E) be positive. One shows that, for all E G J, 

(1.4) & (E) < and > 0. 

To the loop jv,u, we associate the vertical action integral 
S v v defined as 
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Figure 4: The loops for the phases 



(1.5) 



S v ,v{E) 
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where k is a branch of the complex momentum that is continuous on The vertical tunneling 

coefficient is defined to be 



(1.6) 



~S V , V {E) 



The function E 1— > S V)U (E) is real analytic and does not vanish on J. The loop 7^ is oriented so that 
S V>U (E) be positive. 

The index v being chosen as above, we define horizontal action integral u by 



(1.7) 



s h A E ) 



k(C)c*C, 



-/ft. 



where k is a branch of the complex momentum that is continuous on jhv The function E 1— > Sh^ u {E) 
is real analytic and does not vanish on J. The loop 7/^ is oriented so that Sh^ u {E) be positive. The 
horizontal tunneling coefficient is defined as 



(1. 



t hiU (E)=exv[--S h ,„(E) 



In section 10.4, we check that 

(1-9) S h>0 (E) = S h , v (E) and t hj0 (E) = t h>7r (E) . 

One defines 

(1.10) S h (E) = S h , (E) + S h)7r (E) and t h {E) = t hfi (E) ■ t hj7T (E). 

In (1.3), (1.5), and (1.7), only the sign of the integral depends on the choice of the branch of k; this 
sign was fixed by orienting the integration contour; for more details, see sections 10.1 and 10.2. 

1.4. Ergodic family. Before discussing the spectral properties of H ZjE , we recall some general well 
known results from the spectral theory of ergodic operators. 

As 2ir/e is supposed to be irrational, the function x i— ► V(x — z) + a cos(ex) is quasi-periodic in x, and 
the operators defined by (0.1) form an ergodic family (see [22]). 

The ergodicity implies the following consequences: 

(1) the spectrum of H Z;E is almost surely independent of z ([1, 23]); 

(2) the absolutely continuous spectrum and the singular spectrum are almost surely independent 
of z ([23, 18]); 

(3) the discrete spectrum is empty ([23]); 

(4) the Lyapunov exponent exists for almost all z and is independent of z ([23]); it is defined in 
the following way: let x i— ► ip(x) be the solution to the Cauchy problem 

H t , E tp = Eip, W=o = 0, ^'\x=o = 1) 
the following limit (when it exists) defines the Lyapunov exponent: 



log(y/\i/>(x,E,z)\* + W(x,E,z)\i 
Q(E) = @(E,e) := lim * — 



(5) the absolutely continuous spectrum is the essential closure of the set of E where @(E) = 
(the Ishii-Pastur-Kotani Theorem, see [23]); 

(6) the density of states exists for almost all z and is independent of z ([23]); it is defined in the 
following way: for L > 0, let H ZjE .l be the operator H ZjE restricted to the interval [-L, L] with 
the Dirichlet boundary conditions; for E £ R; the following limit (when it exists) defines the 
density of states: 

Ar/J _,> Ar ,^ . .. #{ eigenvalues of H Z£ . L \ess then or equal toE} 
N(E) = N(E,e) := hm — ; 

L— >+oo ZLi 

(7) the density of states is non decreasing; the spectrum of H Z £ is the set of points of increase of 
the density of states. 

1.5. A coarse description of the location of the spectrum in J. Henceforth, we assume that the 
assumptions (H) and (O) are satisfied and that J is a compact interval satisfying (TIBM). Moreover, 
we suppose that 

(T): 2ir ■ minmin(ImC2n-2(-E'), ImC2n+3(^)) > maxmax(5/ l (£'), S Vt o(E), S V ^(E)). 

Note that (T) is verified if the spectrum of Hq has two successive bands that are sufficiently close 
to each other and sufficiently far away from the remainder of the spectrum (this can be checked 
numerically on simple examples, see section 1.8). In section 1.9, we will discuss this assumption 
further. 

Define 

(1.11) 5 := \ mimm(S h (E),S v , (E),S v , w (E)) > 0. 
We prove 

Theorem 1.1. Fix E* £ J. For e sufficiently small, there exists V* C C, a neighborhood of E*, and 
two real analytic functions E i— ► <&o(E,e) and E i— > $ n (E,e), defined on V* satisfying the uniform 
asymptotics 

(1.12) $ (E,e)=$o(E)+o(e), & v (E,e) = K(E) + o(e) where sup |e _1 o(e)| 0, 

Eev* 
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then, for all z, the spectrum of H z e in J Pi V* is contained in the union of the intervals 



jW : = 4° + [-e- 5 °/ £ , e~ s °/ £ ] and if) := E P + [-e~ 5o/£ , e^ ^] 



that is 



^, £ )nJnKc(u4°)ufU4° 



In the sequel, to alleviate the notations, we omit the reference to e in the functions <I>o an d ^tt- 

By (1.4) and (1.12), there exists C > such that, for s sufficiently small, the points defined in (1.13) 
satisfy 



Moreover, for v G {0,7r}, in the interval Jfl V*, the number of points E„ is of order 1/e. 

In the sequel, we refer to the points Eq (resp. i4 ), and, by extension, to the intervals Iq (resp. 
4^) attached to them, as of type (resp. type tt). 

By (1.14) and (1.15), the intervals of type (resp. tt) are two by two disjoints; any interval of type 
(resp. tt) intersects at most a single interval of type tt (resp. 0). 

1.6. A precise description of the location of the spectrum in J. We now describe the spectrum 
of H z e in the intervals defined in Theorem 1.1. Let us assume the interval under consideration is of 
type tt. One needs to distinguish two cases whether this interval intersects or not an interval of type 
0. The intervals of one of the families that do not intersect any interval of the other family are called 
non-resonant, the others being the resonant intervals. 

In the present paper, we only study the non-resonant intervals; the resonant one are studied in [7]. 
The non-resonant is the simplest of the two cases; nevertheless, one already sees that new spectral 
phenomena occur. 

Remark 1.2. One may wonder whether resonances occur. They do occur. Recall that the derivatives 
of $ n and $0 are of opposite signs on J, see (1.4). Hence, as e decreases, on J, the points of type 
and tt move toward each other (at least, in the leading order in e). The motion being continuous, 
they meet. 

Nevertheless, for a generic V, there are only a few resonant intervals in J. On the other hand, for 

symmetric V , there may be numerous resonant energies; e.g., if V is even, then the sequences (E^)i 
(i f ) 

and (E^ )// coincide and all the intervals are resonant! This is due to the fact that the cosine is even; 
it is not true if acos(-) is replaced by a generic potential. 

We will describe our results for the intervals of type tt; mutandi mutandis, the results for the intervals 
of type are the same. One has 

Theorem 1.2. Assume the conditions of Theorem 1.1 are satisfied. Fore sufficiently small, let (Iq )i> 
and (4 )z be the finite sequences of intervals defined in Theorem 1.1. Consider I such that, for any I', 
4 nig = 0* Then, the spectrum of H Zj£ in 4 is contained 4 ? the interval centered at the point 



(1.15) 



(1.14) 



i e < 4° " ^ < Ce, 
±e<E®-E%-*<Ce. 



and of length 
(1.17) 
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2c 



t h (E. 
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V' 



COS 



(1 + 0(1)). 



/ 



The factor A n (V) is positive, and depends only on V and on n (see section 6.2.1). 

In (1.17), o(l) tends to when e tends to 0, uniformly in E £ and / swc/i £/ia£, /or any I', 



(»') 





The fact that each of the intervals 1% does contain some spectrum follows from 

Theorem 1.3. Let dN e {E) denote the density of states measure of E ze . In the case of Theorem 1.2, 
for any I, one has 



I dNJE) = — . 
J & 2tt 



"Level repulsion". Let Eq be the point in the sequence (Eq^)ii closest to E w := E%' . Analyzing 
formulae (1.16) and (1.17), one notices a repulsion between the intervals Iq and I n . 
Indeed, consider the second term in the right hand side of (1.16). As & n (E) > 0, this term has 

the same sign as tan ( ®°( Ew ) ) , Assume that Eq and E v are sufficiently close to each other. As, by 



,(0 



1 - 7T 

definition, — &q(Eq) = — mod tt and as &q(E) < 0, the second term in the right hand side of (1.16) 

is negative (resp. positive) if E n is to the left (resp. right) of Eq. So, there is a repulsion between Iq 
and As the distance from E n to Eq controls the factor 

cos 



the smaller this distance, the larger the repulsion. 

1.7. The Lyapunov exponent and the nature of the spectrum in J. Here, we discuss the 
nature of the spectrum in the interval ij} ■ Therefore, we define 



(1.18) 



A, 



where, for a set A, dist (E, A) denotes the Euclidean distance from E to A. 
1.7.1. The Lyapunov exponent. We prove 



(1.19) 



f(0 

e(E,e) 



Theorem 1.4. On the interval 1^ , the Lyapunov exponent has the following asymptotic 

e 



2tt 



log+ X n (E^)+o(l) 



where o(l) tends to when e tends to 0, uniformly in E € 1$ and I such that, for any i',i?mP = 0. 
Here, log + = max(0,log). 

1.7.2. Sharp drops of the Lyapunov exponent due to the resonance interaction. Formula (1.19) shows 
that the Lyapunov exponent becomes "abnormally small" on the interval I^P when it becomes close 
to one of the points {E Q }. Let us discuss this in more details. 



Assume that (Sh — S Vj7r )(E^p) > 0. If dist {^E^ 1 , [J^Eq '}j > e N (where N is a fixed positive integer) 
then, Theorem 1.4 and formula (1.18) imply that 

&(E, e) = ^-(S h - S V , V ){E®) + when e - 0. 
On the other hand, when E$ is only at a distance of size e _<5//£ (for < 5 < (Sh — S V)7T ) + ) from the 



set of energies {Eq on 1%' , one has 

l r, 



f(0 



@(E,s) = — [(S h -S v ,*)(E®)-6 



+ o(l) when s — > 0. 



Hence, the value of G on ffl drops sharply when E$ approaches the sequence (Eq 

1.7.3. Singular spectrum. As a natural consequence of Theorem 1.4 and the Ishii-Pastur-Kotani The- 
orem [23], we obtain the 

Corollary 1.1. Fix c > 0. For e sufficiently small, if 1% is non-resonant and if e log \ W (E^') > c, 
then, the interval l!p defined in Theorem 1.2 only contains singular spectrum. 

1.7 .4. Absolutely continuous spectrum. If is small on the interval /jj- , most of this interval is made 
of absolutely continuous spectrum; one shows 

Theorem 1.5. For c > 0, there exists rj, a positive constant, and a set of Diophantine numbers 
D C (0, 1) such that 

• asymptotically, D has total measure i.e. 

(1.20) mes(D £ n(0 ' £)) = l + e-^o(l). 

• for e € D sufficiently small, if 1$ is non-resonant and if e log X 7T (E^) < —c, then, one has 

(1.21) =l + o(l), 

mes (Iw ) 

and S ac denotes the absolutely continuous spectrum of H Z £ . 

In (1.20) and (1.21), o(l) tends to when e tends to 0, uniformly in E £ 1$ and I such that, 
for any I', I® n if ) = 0. 

1.7.5. A remark. The nature of the spectrum depends on the interplay between the values of the 
actions S^, S Vt o, S V)n . So, when analyzing our results, it is helpful to keep in mind the following 
observation. As underlined at the end of section 1.5, choosing e carefully, one can arrange that the 
distance between the sequences of energies of type and it be arbitrarily small; moreover, this can be 
done in any compact subinterval of J of length at least Ce (if C is chosen sufficiently large). On such 
an interval, the actions E \— > Sh(E), E >— > S Vt o(E) and E \— > S Vj7T (E) vary at most of Ce. Hence, at the 
expense of choosing e sufficiently small in the right way, we may essentially suppose that there exists 
an energy of type and one of type ir at an arbitrarily small distance from each other such that, on 
an e-neighborhood of these points, the triple E \— > (Sh(E), S Vt o(E), S Vtn (E)) takes any of its possible 

values on J. This means that one can pick the values of eP - E^ and (S h (E), S V)0 (E), S V ^(E)) 
essentially independently of each other. 

Now, let us discuss two new spectral phenomena that can occur under the hypothesis (TIBM). 

1.7.6. Transitions due to the proximity to a resonance. The nature of the spectrum on the intervals 
defined in Theorem 1.2 depends on their distance to the intervals of the other family. The interaction 
can be strong enough to actually change the nature of the spectrum. Let us consider a simple example. 
Assume the interval J satisfies: 

(1.22) mmSh(E) > max max S v v (E) , 

E£j v£{0,tt} e&j ' 

and 

3 

(1.23) - min min S v U (E) > max Sh(E). 

2 u£{0,tt}E£J ' E£j 

Condition (1.22) guarantees that 5q = — min minS v V (E). Hence, there exists c > such that, for 

2 vg{o,7t} 

E € J and v £ {0, 7r}, 

(1.24) S h (E) - S V>V {E) - S < -c < 0. 

Consider now Ip and l!p both non resonant located in Jfl 14. Then, 

• if the two intervals are distant of at least e N (where iV is a fixed integer) from each other, 
condition (1.22) guarantees that, on these intervals, the spectrum is controlled by Corollary 1.1 
and its analogue for the intervals of type 0. 
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• if the two intervals are distant of Cexp (— 5q/e) from each other, then condition (1.24) guar- 
antees that, on these intervals, the spectrum is controlled by Theorem 1.5 and its analogue for 
the intervals of type 0. 

That intervals J where both (1.22) and (1.23) hold exist can be checked numerically, see section 1.8. 
Thus, not only does the location of the spectrum depend of the distance separating intervals of type 
for neighboring intervals of type ir, but so does also the nature of the spectrum. Transition can occur 
due to this interaction phenomenon: spectrum that would be singular were the intervals sufficiently 
distant from each other can become absolutely continuous when they are close to each other (see 
Fig. 5(a)). 

1.7.7. Alternating spectra. To describe this phenomenon, to keep things simple, assume that, in V*H J, 
the distance between the points {£q^} and the points {E^ ^} is larger than e N (for some fixed N); 
hence, all energies are non-resonant in V* H J. Taking Theorem 1.5 and Corollary 1.1 into account, 
we see that, on /q (resp. 1$ the nature of the spectrum is determined by the size of the ratio 
t Vt0 {E$ ] )/t h (E®) (resp. t v ^(Ep)/t h (E { J' ] )). So, if for some 5 > 0, one has 

(1.25) V£ g J n K, S ViV (E) - S h {E) >8 et S vfi (E) - S h (E) < -5, 

then, in V* H J, the sequences of type and ir contain spectra of "opposite" nature: the spectrum in 
the intervals of type is singular, and that in the intervals of type ir is, mostly, absolutely continuous. 
This holds under the Diophantine condition on e spelled out in Theorem 1.5. Hence, one obtains an 
interlacing of intervals containing spectra of "opposite" types, see Fig. 5(b). In this case, the number 
of Anderson transitions in V* H J is of order 1/e. 

One can check numerically that the condition (1.25) is fulfilled for some energy region V* and some 
values of a (see section 1.8). 



E 



sing ac - a - c - sing a.c. sing a.c. sing 

(a) Almost resonant transitions (b) Alternating spectra 



Figure 5: Two new spectral phenomena 



1.8. Numerical computations. We now turn to numerical results showing that the multiple phe- 
nomena described in sections 1.7.6 and 1.7.7 do occur. 

All these phenomena only depend on the values of the actions Sh, S Vj o, S V7r . For special potentials 
V, they are quite easy to compute numerically. 

We pick V to be a two-gap potential; for these potentials, the Bloch quasi-momentum k (see sec- 
tion 1.1.2) is explicitly given by a hyper-elliptic integral ([17, 20]). The actions then become easily 
computable. As the spectrum of Hq = —A + V only has two gaps, we write a (Ho) = [E±,E2\ U 
[£3, £4] U [£5, +00 [. In the computations, we take the values 

Ei = 0, £ 2 = 3.8571429, £3 = 6.8571429, £4 = 12.100395, and £5 = 100.70923. 

On the figure 6, we represented the part of the (a, £)-plane where the condition (TIBM) is satisfied 
for n = 1. Its boundary consists of the straight lines £ = E\ + a, £ = £2 + a, £ = £3 — a and 
£ = £4 — a. Denote it by A. 

The computations show that (T) is satisfied in the whole of A. As n = 1, one has £2n-2 = —00. So, 
it suffices to check (T) for C,2n+3 = Cs- (T) can then be understood as a consequence of the fact that 
£5 — £4 is large. 

On the figure 6, one sees that, for non-resonant intervals, 

• the zones where one has alternating spectral types (see section 1.7.7) are those where either 

S v ,o < Sh < S^tt or S v<n < Sh < S V} o 
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E = E± — a 




Figure 6: Comparing the actions 



• the transitions due to the proximity of the resonant situation (see section 1.7.6) take place in 
the part of region {Sh > max(5 Uj7r , S^o)} sufficiently close to {Sh < mm(S V)7T , S Vt o)}. 

1.9. Comments, generalizations and remarks. About assumption (T), its purpose is to select 
which tunneling coefficients play the main role in the spectral behavior of H ze in the interval J; 
this assumptions guarantees that it is the tunneling coefficients associated to the loops defined in 
section 1.3.2 that give rise to the principal terms in the asymptotics of the monodromy matrix that 
we describe in section 2. 

In the present paper, we restricted ourselves to perturbations of Hq of the form a cos. As will be 
seen from the proofs, this is not necessary. The essential special features of the cosine that were used 
are the simplicity of its reciprocal function (that is multivalued on C). More precisely, the assumption 
that is really needed is that the geometry of the objects of the complex WKB method that is used to 
compute the asymptotics of the monodromy matrix be as simple as that for the cosine. This geometry 
does not only depend on the perturbation; it also depends on the interval J under consideration and 
on the Bloch quasi- momentum of Hq. The precise assumptions needed to have our analysis work are 
requirements on the conformal properties of the complex momentum. 

The methods developed in [11, 12, 8, 9, 10] are quite general; using them, one can certainly 
analyze more complicated situations i.e., more general adiabatic perturbations of Hq. Nevertheless, 
the computations may become much more complicated than those found in the present paper. 

1.10. Asymptotic notations. We now define some notations that will be used throughout the paper. 
Below C denotes different positive constants independent of e, E and E n . 

When writing / = 0(g), we mean that there exits C > such that |/| < C\g\ for all e, z, E in 
consideration. 

When writing / = o(g), we mean that there exists e i— > c(e), a function such that 

• |/| < c ( e )l<?l f° r an e > z -> E m consideration; 

• c(e) — ► when e — > 0. 

When writing / x g, we mean that there exists C > 1 such that C -1 |g| < |/| < C\g\ for all e, z, E in 
consideration. 

When writing error estimates, the symbol 0(/i, f2, ■ ■ ■ fn) denotes functions satisfying the estimate 

(1-26) \0(f u f 2 ,... f n )\ < C(\h\ + \f 2 \ + .-. |/n|), 

with a positive constant C independent of z, E and e under consideration. 

li 



2. The monodromy matrix 
In this section, we consider the quasi-periodic differential equation 
d 2 

(2.1) — —rijj(x) + (V(x-z) + acos(ex))ip(x)=E?p(x), x G R, 

and recall the definition of the monodromy matrix and of the monodromy equation for (2.1). We 
also recall how these objects are related to the spectral theory of the operator H g>£ defined in (0.1). 
Finally, we describe two monodromy matrices for (2.1). 

2.1. The monodromy matrices and the monodromy equation. We now follow [11, 12], where 
the reader can find more details, results and their proofs. 

2.1.1. The definition of the monodromy matrix. For any z fixed, let (ipj(x, z)) je | 12 j be two linearly 
independent solutions of equation (2.1). We say that they form a consistent basis if their Wronskian 
is independent of z, and, if for j G {1, 2} and all x and z, 

(2.2) ipj(x, z + l) =ipj{x, z). 

As (ijjj(x, z))jg{ lj2 } are solutions to equation (2.1), so are the functions ((x,z) i— > ipj(x + 2ir/e,z + 
27r/e))j e r 12 |. Therefore, one can write 



(2.3) * (x + 2ir/s, z + 2-ir/e) = M (z, J5)*(x, z), V(x,z) 



ipi(x, z) 

i>2{x,z) 



where M (z,E) is a 2 x 2 matrix with coefficients independent of x. The matrix M(z, E) is called the 
monodromy matrix corresponding to the basis {ip j) j&f 1,2}- To simplify the notations, we often drop 
the E dependence when not useful. 

For any consistent basis, the monodromy matrix satisfies 

(2.4) detAT(z) = l, M(z + l) = M(z), \fz. 

2.1.2. The monodromy equation and the link with the spectral theory of H z ^ e . Set 

27T 

(2.5) ft, = — modi. 

Let M be the monodromy matrix corresponding to the consistent basis (^j)j=i,2- Consider the mon- 
odromy equation 

(2.6) F(n+1) = M (z + nh)F(n), where F{n) G C 2 , Vn G Z. 

The spectral properties of H Z £ defined in (0.1) are tightly related to the behavior of solutions of (2.6). 
For now we will give a simple example of this relation; more examples will be given in the course of 
the paper. 

Recall the definition of the Lyapunov exponent for a matrix cocycle. Let z \— > M(z) be an SL(C, 2)- 
valued 1-periodic function of the real variable z. Let h be a positive irrational number. The Lyapunov 
exponent for the matrix cocycle (M, h) is the the limit (when it exists) 

(2.7) 6(M,h)= lim - log \\M(z + Lh) ■ M(z + (L - l)h) ■ ■ ■ M(z + h) ■ M(z)\\. 

L— >+oo L 

Actually, if M is sufficiently regular in z (say, if it belongs to L°°), then 6(M,h) exists for almost 
every z and does not depend on z, see e.g. [23]. 

One has 

Theorem 2.1 ([1 I]). Let h be defined by (2.5). Let z \— > M(z,E) be a monodromy matrix for 
equation (2.1) corresponding to basis solutions that are locally bounded in (x,z) together with their 
derivatives in x. 

The Lyapunov exponents Q(E,e) and 6(M(-,E),h) satisfy the relation 

(2.8) @(E,h) = ^e(M(-,E),h). 
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2.2. Asymptotics of two monodromy matrices. Recall that, in the interval J, the spectrum 
of H z e is contained in two sequences of subintervals of J, see Theorem 1.1. So, we consider two 
monodromy matrices, one for each sequence; for v G {0,7r}, the monodromy matrix M„ is used to 

study the spectrum located near the points (Eiy)i- 

In this section, we first describe the monodromy matrix M n in detail. Then, we briefly discuss the 
monodromy matrix Mq. 

Fix v G {0, 7r}. The monodromy matrix M v is analytic in z and E and has the following structure: 
(2.9) M v : 

where, for (z\, ■ ■ ■ , z n ) \— > g(z\, ■ ■ ■ , z n ), an analytic function, we have defined 



(2.10) g*{zi, ■ • ■ , z n ) = g(zi, ■■■ ,z n ). 

When describing the asymptotics of the monodromy matrices, we use the following notations: 

• for Y > 0, we let 

(2.11) T Y = e- 2nY / £ ; 

• we put 

(2.12) p{z) = e 2 -l Im < 
One has 

Theorem 2.2. There exists V*, a complex neighborhood of E*, such that, for sufficiently small e, the 
following holds. Let 

(2.13) Y m = - 1 - inf m&x(S V)0 (E),S V)7T (E)), Y M = ^- sup m&x(S vfi (E), S V 7r (E), S h (E)). 

ztt .EeJnv; 2ir E^jnv* 

There exists Y > Ym and a consistent basis of solutions of (2.1) for which the monodromy matrix 
(z,E) i— > M n (z,E) is analytic in the domain {z£C: |Imz| < -j} x V* an d has the form (2.9). 
Fix < y < Y m . Let V£ = {E G V* : \ImE\ < e}. In the domain 



(2.14) {zeC: |Imz| < ^} 

the coefficients of M v admit the asymptotic representations: 



x v; 



* ! 



(2.15) A 7T = 2 — h-e e l- + 0\ + o\T h ,—^ — , T vfi p(z), T V)7r p(z ( 

and 

/01 m n a^-e^Co 1 / 1 i*a . a _i*o\ . n T Y p(z) . . , 

(2.16) 5^ = 2 — h 2 g£ \0 e £ " l+Oir fe , — , T vfi p(z), T Vi7T p(z / 

with 

(2.17) C = ^ (aoe^+c^e-- 

In i/iese formulae, for v G {0, 7r} ; (z,E) \— > a u (z,E) is an analytic function and is 1-periodic in z; it 
admits the asymptotics 

(2.18) a v = 1 + T VtU e 2 ™( z - z "( E » + O (T Y p{z)) . 

The quantities E \— > <& V (E), E \— > T V)U (E), E i— ► T^E), E \— > 9(E) and E \— > z v [E) are real analytic 
functions; they are independent of z; for E G V£, they admit the asymptotics: 

(2.19) & v {E)=$ v (E)+o{e), 

(2.20) T h (E) = t h (E)(l + o(l)), T VtV (E) = t v , u (E) (1 + o(l)), 

where Q v andth, t V)U are the phase integrals and the tunneling coefficients defined in section 1.3; 

• 

(2.21) 6(E)=6 n (V){l + o{l)), 

where n (V) is the constant defined in section 6.2; it is positive and depends only on n and V ; 
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(2.22) „ W -» W - *-ro- *°<*> -i -Ml). 

27T£ £ 

(2.23) z' v {E) = 0{l). 

Note that the terms containing 9 in the asymptotics (2.15) and (2.16) are bounded independently of 
e. So, with exponentially high accuracy, the coefficients A n and B n are proportional. 

Remark 2.1. The description of the monodromy matrix Mo is similar to that of M n : in Theorem 2.2, 
one has to change 

(1) the indexes and tt by respectively tt and 0; 

(2) the quantity 9 by 1/0; 

(3) z (E) by z (£) + /t in formulae (2.18). 

Most of the analysis used to construct Mq is the same as that for M n . The differences are described 
in section 5. 

Theorem 2.2 is the central technical result of the paper. In the next two sections, we use Theorem 2.2 
to study the spectrum of H z ^ e , and the remainder of the paper is devoted to its proof. 

2.2.1. Useful observations. We now turn to a collection of estimates used when deriving the results of 
sections 1.5, 1.6 and 1.7 from Theorem 2.2. We begin with 

Lemma 2.1. Let J* C R be a compact interval inside V*. There exists a neighborhood of J*, say V*, 
and C > such that, for sufficiently small e, for E 6 V* and v E {0, tt}, one has 

(2.24) \& v (E)\ + \ty(E)\<C, 
and 

(2-25) ^<\K(E)\- 

Proof. Recall that the phase integrals § u are independent of e, analytic in a neighborhood of J, 
and, on J, the derivatives & V (E) are bounded away from zero, see (1.4). Therefore, the statements of 
Lemma 2.1 follow from (2.19) and the Cauchy estimates for the derivatives of analytic functions (o(e) 
in (2.19) is analytic in the domain V*, and, therefore, on any its fixed compact, one has the uniform 
estimates: ^o(e) = o(e) and ^j^o(e) = o(e)). This completes the proof of Lemma 2.1. □ 
We also prove 

Lemma 2.2. For sufficiently small e, for u £ {0, tt}, in the domain (2.14), one has 

(2.26) a v = l + 0(T VtV p{z)) = l + o(l), 

(2.27) p(z)\T ViU (E)\ = o(l), 

(2.28) e «M*0A x l, 

(2.29) \T h (E)\ + \T V , U (E)\ + T Y < Ce~ 2& ^\ 

(2.30) T Y = o(T h (E)) and Ty = o(T VyU (E)), 

(2.31) Ce- 2nYM/£ < \T h {E)\ and I e - 2 -W £ < \T V>U (E)\ < Ce- 2wYm/£ , 

(2.32) \9{E)\ x 1, 

(2.33) |e 27ri ^|xl. 
All the above estimates are uniform. 

Proof. As z v is real analytic, estimate (2.33) follows from (2.23) and the definition of V £ . Esti- 
mate (2.32) follows from (2.21) as 9 n (V) is a positive constant depending only on n and V. The 
estimates (2.31) follow from (2.20) and the definitions of the tunneling coefficients, of the domain V^f 
and numbers Y m and Ym- Estimates (2.30) follow from (2.31) as Y > Ym- Estimates (2.29) follow 
from (2.30), (2.20) and the definition of <5o- Estimate (2.27) follows from (2.31) as in the domain (2.14), 
one has |Imz| < y/e, and y < Y m . Estimate (2.28) follows from (2.24), the definition of the domain 
V£ and from the real analyticity of the phase integrals. The inequalities in (2.26) follow from (2.30) 
and (2.27). This completes the proof of Lemma 2.2. □ 
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3. Rough characterization of the spectrum and a new monodromy matrix 

In this section, we first obtain a rough description of the location of the spectrum of H Z)£ i.e., we 
prove Theorem 1.1. Then, we change the consistent basis so that, in a neighborhood of the spectrum, 
the new monodromy matrix have a form more convenient for the spectral study. 

3.1. The scalar equation. Our analysis of the spectrum is based on the analysis of solutions of the 
monodromy equation with the monodromy matrices described in the previous section. A monodromy 
equation is a first order finite difference 2-dimensional system of equations, see (2.3). Instead, of 
working with this system, we study an equivalent scalar second order finite difference equation. To 
derive this equation, we use the following elementary observation 

Lemma 3.1. Let M : z *—> M(z) be a SL(2,C)-valued matrix function of the real variable z, and let 
h be a real number. Assume that M\2{z) ^ for all z. Define 

(3.1) p(z) = M 12 (z)/M 12 (z - h), v(z) = M u (z) + p(z) M 22 (z). 

A fucntion \Pi : Z — > C is the first component of a vector function \E' : Z — > C 2 satisfying the equation 

*(fc + 1) = M(hk + z)t>(k), V/cGZ, 
if and only if it satisfies the equation 

(3.2) Vi(k + l) + p(hk + z)Vi(k-l) =v(hk + z)V 1 {k), V/c G Z. 

The reduction from the monodromy equation to the scalar equations (3.2) has already been used in [4] 
and [12]. To characterize the location of the spectrum of (0.1), we use 

Proposition 3.1. Fix E in equation (2.1). Let f and g form a consistent basis in the space of the 
solutions of (2.1), and let M be the corresponding monodromy matrix. 

Assume that the functions (x, z) \— > f(x, z), (x, z) \— > g(x, z), (x, z) \— > d x f(x, z) and (x, z) \— > d x g(x, z) 
are continuous on M 2 . 

Suppose that min|Mi2(z)l > 0. In terms of M, define the functions p and v by (3.1) and define h 
by (2.5). Let 

( 1 \ 2 

(3.3) max \p(z)\ < - min \v(z)\ , indp= indv = 0, 

zeR \2 zeR J 

where indg is the index of a continuous periodic function g. 
Then, E is in the resolvent set of (0.1). 

The proof of this proposition immediately follows from Proposition 4.1 and Lemma 4.1 in [12] based 
on the analysis in [4]. 

Remark 3.1. This proposition is very effective if the coefficient M\ 2 of the monodromy matrix is 
close to a constant. Then, it roughly says that the spectrum is located in the intervals where the 
absolute value of the trace of the monodromy matrix is larger than 2. This is the condition one meets 
in the classical theory of the periodic Schrodinger operator ([6]). 

3.2. Rough characterization of the location of the spectrum. We now prove Theorem 1.1. 
Pick E* G J. Let V* be as in Theorem 2.2. Consider the sequences (Ei l) )i and (E$\ defined by the 
quantization conditions (1.13). 

Introduce 5o by (1.11). Let J* be a compact subinterval of J Pi V*. One has 

Lemma 3.2. Pick < a < 1. For e sufficiently small, in J*, the spectrum of H Z)£ is contained in the 
e a e~ s °^ £ -neighborhood of the points (E$)i and (E^) m defined by the quantization conditions (1.13). 

Lemma 3.2 implies Theorem 1.1 at the possible expense of reducing 14 somewhat. 
Proof. Define 

(3.4) y rough = {E G J* : \E- E^ m) \ > e a e~ s ° /£ , Vm}. 

We shall prove that, for e small enough, the spectrum of H Z)£ in T4 oug h is contained in the e a e~ s ^ e - 
neighborhood of the points (eI^)i. 

In the remainder of this proof, we assume that e is sufficiently small for the statements of Theorem 2.2 
and Lemma 2.1 to hold. 
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The proof consists of the following steps. 
1. We prove that, for e sufficiently small, 



inf 



cos 



ME) 



> e 



So/e 



This follows from (3.4), from the definition of the set {Eq }, and from (2.25). 
2. We check that, for E G J*, and for z£l, each of the functions A v and B n has the form 

2 



e^ /£ cos(4> /e) + 0(e 



■2S /e\ 



Indeed, by the first inequality from (2.26), for v £ {0,7r}, z£l and E £ J*, one has 

ol v = l + 0(T VjU ). 

By means of this estimate and of (2.28) and (2.32), we transform the right hand sides both in (2.15) 
and (2.16) to the form 

2 



j*«' e cos(4> /e) + 0(T vfil T„ j7r , T h , T Y ) 



This and (2.29) imply that A v and B n have the requested form. 

3. Let (z,E) i — > p(z,E) be the function defined by (3.1) for M = M n (z,E). The previous two steps 
imply that there exists C > such that, for e sufficiently small, one has 

sup sup \p(z,E) - 1| < Ce~ s ° /e . 
zeR EeV roxtgb 

4. Let (z, E) \—* v(z, E) be the function defined by (3.1) for M = M n (z, E). The previous three steps 
imply that, for Q £ M and E G 14ough 5 one has 



u(^£0=Ar + < + (p(«,^0-lK 

2 



2/, 



2cos(^/e) cos($ /e) + 0{e~ 2So/£ ) 



+ 



e^ /e cos($ /e) + 0(e- 2So/e )) 0(e 



-<5o/ £N 



cos($ /e) (cos($ n /e) + 0(e _,5 ° /e ) 



5. There exists C > such that, for e sufficiently small, if E £ a(H z>£ ) n Plough) then 



(3.5) 



cos 









cos 



goCg) 



Indeed, by steps 1 and 2, for sufficiently small e, for £7 G Koughi one has 

min|J3 w (z,£0| > 0, vadBJ-, E) = 0. 

Moreover, by steps 3 and 4, there exists C > such that, for e sufficiently small, for E £ T4 oug h, if 



cos 



> C e~ s °/ £ + 



■i-o(g) / ' 



then, one has 



min \v(z, E)\ 2 > 4max \p(z, E)\, indf (•, E) = 0. 



These two observations and Proposition 3.1 complete the proof of (3.5). 
6. In view of (2.29) and of the first step, inequality (3.5) implies that 

\cos($ n (E)/e)\ < Ce~ 5o ' £ . 

By the definition of (E% )i and Lemma 2.1, this implies that there exists I such that \E — E%'\ < 
Cee~ s °/ e . This completes the proof of Lemma 3.2. □ 
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3.3. A new monodromy matrix. As said, to study the spectrum, instead of working with the 
monodromy equation itself, it is more convenient to work with the equivalent scalar equation (3.2). 
The use of this equation is very effective when M12, the element of the monodromy matrix, is close to a 
constant, and M\\ (or/and its derivative in E) is much larger than M22- To satisfy these requirements 

for E near the points (E$)i, we introduce a new monodromy matrix. Therefore, we make the following 
simple observation: 

Lemma 3.3. Recall that h is defined by (2.5). Let M be a monodromy matrix for equation (2.1), and 
let U : z 1— > U{z) £ SL(2,C) be a 1-periodic matrix function. Then, 

(3.6) M u (z) = U{z + h)M(z)U{z)~ 1 

is also a monodromy matrix for equation (2.1). 

Proof. Let f\ and fi be the solutions of (2.1) that form a consistent basis for which M is the 
monodromy matrix. The components of the vector 

J2(X,Z) 

are also solutions of (2.1); they form a consistent basis, and M u is the corresponding monodromy 
matrix. □ 



(3.7) 



T{x, z) = U{z)F(x, z), F(x, z) 



For (z,E) in the domain (2.14), we define the new monodromy matrix M u choosing M = M w (z,E), 
the matrix described in Theorem 2.2, and 

1/1 1\ fj(z) 



(3. 



U(z) 



2 \-i ij V 7*(*)/ ' K-T-v \/ a^z) 

Recall that, for (z, E) being in the domain (2.14), by Lemma 2.2, one has a = 1 + o(l) when e tends 
to 0. So, we define a branch of 7 analytic in this domain by the condition = 1 + o(l). 

Then, one proves 



Theorem 3.1. In the case of Theorem 2.2, in the domain (2.14), the monodromy matrix (z,E) 
M u (z,E) is real analytic and admits the representation: 



(3.9) 

where 

(3.10) 



M u (z, E) = P(z, E) + Q(z, E)+0[ T h , p(z)^-, p(z)T vfi , p(z)T V)V ) , 



Ty 



P(z,E) 



4 fC n (z,E)C (z,E) -S n (z,E)C (z,E) 











(3.11) Q(z,E) 
In these formulae 



1 C os £*=to + e cos ^ cos ^ 

o e s e 

$0 



-0 sin ^ 



l sint^to _#sin^ cos to 

f £ „ £ £ 

6* sin ^ sin ^ 



(3.12) 



Ctt 



1 

2i 



a- e 



-i*7r/£ 



and (2, .E) i — > d(z, .E) is an analytic function that admits the asymptotics: 

(3.13) a* = 1 + T„ )7r [cos(2tt(z - z„)) + isin(27r(z - /» - z^))] + 0(p 2 (z)Tl w , p{z)T Y ). 

All the above estimates are uniform in the domain (2.14). 

Proof. The monodromy matrix M u is analytic in the domain (2.14) as M n and U are. As the 
consistent basis in Theorem 2.2 consists of a pair of solutions of the form f\ = f and fi = /*, for 
U given by (3.8), formula (3.7) defines two consistent solutions of (2.1), say fY and fif, such that, 
for x fixed, (z,E) 1— ► /f 7 (x, z, E) and (z,E) i— > (x, z, E) are real analytic. So, the new monodromy 
matrix (z,E) \—> M u (z,E) is also real analytic. 
Compute Mg. By (3.8) and (3.6), 

5 + 5* 



(3.14) 



where 5 = 7 *(z) [ 7 (z + /i)^(z) + 7 *(z + . 
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The definition of 7 yields 



s _ l(z + h) 



7(2) 

Substituting the asymptotic representations (2.15) and (2.16) into this expression, and using the real 
analyticity of Th, $0) ^Vj 9 and Co, we get 

(3.15) _l ^Z[M ^ + fe ) f!!^! f I e -i*o/e , e i*o/e 

7 (z) 2 
7(z + /t) c c ai$,/ e 7(g + fe) a^(z) 
7(z) 7(2;) a*(z) 

where d,r(z) = 7 ^ +fe ) ^(-z), and denotes OCT^, Typiz) /Ty n T v op(z), T V7T p(z)). By the estimates 
of Lemma 2.2, from (3.15), one obtains 

S = ^- a n {z) j**/ e C (z) + e^/ £ + eos($ /e)) + O. 

Substituting this result into (3.14), we get the formula announced for M|{ in Theorem 3.1. 

The other coefficients of the matrix M u are computed analogously; so, we omit the details. 

To complete the proof of Theorem 3.1, it remains only to check (3.13). Put a n< i = a n — I. By 

Lemma 2.2, one has a W) i = 0(pT V)7T ). Therefore, 

(3.16) 1{Z) V <{z-h) ) 

= 1 + 2 + + «->i(^ - ^ - " ^)) + 0{p 2 (z)Tl^). 

In view of (2.18), one has 0^,1 = T ViW e M (*-**( E fi+0 (T Y p(z)). Substituting this in (3.16) yields (3.13). 
This completes the proof of Theorem 3.1. □ 

Finally, we note that, similarly to (2.26), one proves that 
Lemma 3.4. Uniformly in (z,E) in the domain (2.14), one has 

a„ = l + 0(T VtV p(z)) = 1 + 0{e~ 2 ^ Y ™-yV e ) = 1 + o(l). 

4. The spectrum in the "non-resonant" case 

We now prove the results on the spectrum of H Zj£ formulated in Theorems 1.2, 1.3, 1.4, 1.5 and 
Corollary 1.1. 

Pick E* G J. Let V* be as in Theorem 2.2. Let J* C V* PI R be a compact interval centered at E*. 
We always assume that e is so small that the statements of Theorem 3.1 and Lemma 2.1 hold. 
Let E n be one of the points of (E$)i in J*. We assume that E n satisfy the non resonant condition 

(4.1) inf \E^ m) - EJ > 2e~ 5 » /£ . 

m 

In this section, we fix a satisfying 

< a < 1. 

and study the spectrum in the e a e _<5o//£ -neighborhood of E n . 

Our main tool will be the scalar equation (3.2); recall that we consider the one associated to the 
monodromy matrix M u described in Theorem 3.1. 

In the sequel, we use the notations defined in section 1.10. Now, all the symbols are uniform in 

En- 

4.1. Coefficients of the scalar equation. Here, we analyze the coefficients of the scalar equation 
for energies E satisfying 

(4.2) \E-E n \ < e a e~ 5o/£ . 

18 



4.1.1. The results. We now define a few objects that we shall use throughout the analysis. Let 



(4.3) o~ n = — sin 



As E n G {E^}, one has either a n = +1 or a w = — 1. 
Let 



(4.4) F,(E) = ^ * cos ( MM) VAM {E _ E ,) _ 2K(V) sin < '*«< £ - 



The factor A n (V) is defined in (6.3). The coefficient F n (E) will play the role of an "effective spectral 
parameter" . 

Also, we define the factor 

(Art \ A T "A^) fMEn) 

(4.5) Ajr = 4a n cos 1 



This factor will play the role of an "effective coupling constant" . Finally, we let 



(4.6) £i = min i S , [ 2ttY - max S h (E) 

I \ EeJnv* 

where 5q is defined by (1.11) and Y is the constant from Theorem 2.2. We note that 

<5i Ym 
< — < — . 
2vr - 2 

These inequalities follow from the inequalities Y > Ym and 5o < irY m in which Ym and Y m are the 
numbers defined by (2.13). 
We prove 

Proposition 4.1. Let p u and v u be the the coefficients p and v of the scalar equation (3.2) corre- 
sponding to the monodromy matrix M u . 
Assume that the condition (4.1) is satisfied. 

Fix < y < 5\/{2it). Then, the strip {|Imz| < y/e}, for E satisfying (4.2), one has M\2 7^ 0, and 
the coefficients p u and v u admit the following asymptotic representations 

(4.7) p u (z, E) = l + (p{z)ee- So/£ 

v u (z, E) = \f{E) + K sin(2vr(z -h- z^{E v ))) + o(p(z)X 7r (E)) + O (p(z)e~ 5l / £ ) } 

(4-8) 1 ( ( \i 

• M + O (p{z)ee~ Sa/£ 

Here, the function E i— > F(E) is independent of z; F(E) and F'(E) admit the asymptotic representa- 
tions: 

(4.9) F(E) = F 7T (E)(l + o(l))+o(l), and F'(E) = F^(E)(1 + o(l)). 

We often shall use simplified versions of (4.7) and (4.8), namely 
Corollary 4.1. In the case of Proposition (4.1), one has 

(4.10) p u (z,E) = l + o(l), 

(4.11) v u (z, E) = {F n (E) + cos(27r(z - h - Z7T (E n ))) + o{Kp(z)) + o(l)} (1 + o(l)). 
Proof. For < y < Si/(2n) and |Imz| < y/e, one has 

(4.12) p{z)e- So/£ + p(z)e- Sl/£ < e -( s ^y)/c. 

Representation (4.10) is obtained from (4.7) by means of (4.12). Representation (4.11) is obtained 
from (4.8) by means of (4.12) and (4.9). □ 
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4.1.2. Proof of Proposition 4-1- Let us begin with three lemmas. First, we collect simple observations 
following from Taylor's formula. One has 

Lemma 4.1. For e sufficiently small, for all E satisfying (4.2), for u £ {0,7r}, one has 

(4.13) |cos (^{E)/e) \ < Ce a ~ l e~ s °l £ , 

(4.14) co S ($ n (E)/e) =a 7T e- 1 K(E 7 r)(E-E 7T )[l + 0(e a - 1 e- s °/ e )], 

(4.15) sin ($„(E)/e) = sin ($ v (E n )/e) + O^-V 4 "/"), 

(4.16) |cos ($ (E)/e) \ > Ce^e^ ^ 

(4.17) co^ {E)/e) = cob@ (E v )/ £ ) (1 + 0{e a )) , 

(4.18) T h (E) = T h (E n )(l + 0(e a - l e- 5 °/ £ )), T V , U (E) = T v , v {E„)(l + 0{e a ~ l e~ & »' £ )). 

Proof. These results follow from the Taylor formula. When proving the first five results, one uses (2.24) 
and (2.25) and has to keep in mind the definitions of Eq and E n . We omit the elementary details. 
The two estimates (4.18) are proved in one and the same way. We prove only the first one. Therefore, 
one uses the Taylor formula for log Th (E) in the neighborhood (4.2) of E w . By (2.20) and the definition 
of th, one has logTh(E) = — ^Sh{E) + g(E), where g(E) = o(l) uniformly in V*. The estimates 
ISfJ-E')! < C and ^ = o(l) hold uniformly on any fixed compact of V* (the last estimate follows from 
the Cauchy estimates). This implies that, for E in a fixed compact of V*, 

(4.19) \T' h {E)\ < Ce- l \T h {E% 

and this estimate implies the estimate for Th from (4.18). This completes the proof of Lemma 4.1. □ 
We also prepare simplified representations for factors Co, S n and C n defined in (2.17) and (3.12). We 
prove 

Lemma 4.2. Fix y as in Proposition 4-1- Under condition (4.1), for \lmz\ < y/e and E satisfy- 
ing (4.2), one has 



(4.20) C = cos {ME)/e) (1 + (p(*)«T* /e )), 

(4.21) C w = cos (K{E)/e) + a^{E^) sin(2^(z -h- z^E*)) + oipT^E^), 

(4.22) S n = sin {^(E)/e) (1 + 0(pe~ 2S °/ £ )). 

Proof. The definitions of C and S n , (2.17) and (3.12), and (2.28), (2.26) imply that 

(4.23) C = cos {$ (E)/e) + 0{pT vfi ) and S n = sin ($ n (E)/e) + 0(pT v>7r ). 

Representation (4.20) follows from (4.23), from estimate (4.16) and from (2.29). Similarly, (4.22) 
follows from (4.23), (4.13) and (2.29). 

Prove (4.21). The definition of C n , (3.12), and representation (3.13) imply that 

C n = cos C + T V ^{E) Tcos C h®\ c(z) - sin ( ^1) s(z)] + Q(j?T 2 pTy). 



where s(z) = sin(2-7r(z — h — z n (E)) and c(z) = cos(2tt(z — z n (E)). Now, representation (4.21) follows 
from (4.13), from (4.18) and from estimates (2.30) and (2.27). □ 

Turn to the proof of Proposition 4.1. Compute p . By (3.9), we have 
(4.24) Mf 2 = P 12 + Q 12 + R l2 , R 12 = O [T h , p{z)^, p(z)T vfi , p(z)T V)7T 



Show that, for E satisfying (4.2) and |Imz| < y/e, one has 

(4.25) P 12 = i- sm(^(E)/e) co S (<b (E)/e) (l + 0(p(z)ee- 5 ^ £ )) , |Q 12 | + \R l2 \ < C. 

The estimate for P\ 2 follows from Lemma 4.2. The estimate for Q\ 2 follows from (2.32) and (2.28). 
Check the estimate for R\ 2 . By (2.29) and the definition of <5i, one has |JRi 2 | < Cp(z)e~ Sl / £ . Recall 
that p = e 27r l Imz L As y < 5i/(2tt), for |Imz| < y/e we get 

(4.26) p(z)e- Sl/E < e -( 5 i- 2 ^)A < C. 
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This implies the announced estimate for R\ 2 and completes the proof of (4.25). 

For E satisfying (4.2), as satisfies (1.13), for e sufficiently small, one has | sin($ 7r (E)/e) | > 1/2 ; 
taking (2.29) and (4.16) into account, we get 

1 -1 



— sin(3v(£)/ e )cos($o(£)/e) 



< Cee~ &a/e . 



From this, (4.24) and (4.25), one deduces 

4 



(4.27) 



<2 



n 



sin(^(E)/e) cos{<S> (E)/e) ( 1 + 0{p{z)ee- So/e ) ) . 



In view of (4.26), there exists £q > such that, for < e < £q, the- error term in (4.27) be smaller 
than 1/2. From now on, we assume that < e < £q. Then, we get 7^ 0, and, as p u (z) = 
M^ 2 {z)/M^ 2 {z - h), the representation (4.27) implies (4.7). 

Now, let us compute v u . Note that v u (z, E) = Afg (z, E)+M% 2 (z -h,E) + (p u (z,E)- 1)M% 2 (z -h,E). 
Using the representations (3.9), (3.10) and (3.11), we transform this expression to 

(4.28) v u (z,E)=P u (z,E) + (Q 11 (E)+Q 22 (E))+R(z,E), 

(4.29) R = (p u (z, E) - IXQ22GE) + ri(z, E)) + r 2 (z, E), 

(4.30) rj (z,E) = O (T h , p(z)3fc, p(z)T vfi , p(z)T v ^ for j G {1,2}. 
We now show that 



(4.31) 
where 

(4.32) F(E) 

and that 

(4.33) 
(4.34) 



P u (z,E) = {F(E) + X JT s(z) + o(pX n )) (l + g(z,E)) 



4c os^cos^ 



T h (E) 



s{z) = sin(2^(z - h - z^E*))), \g(z, E)\ < C p(z)ee~ So/e , 



\Qn{E)\ + IQ22OEOI < C, \R(z,E)\ < Cp{z)e~ 5 ^^ 
\g(z,E)\<Ce, \R(z,E)\<C. 



Lemma 4.2 implies that 
(4.35) 



ni{Z,&) jrjE) 



_4_ . cos MS (! + o (p(z)ee- s °/z)) (cos ^ + a n T Vi7V s(z) + o(p(z)T v 
F{E) + \{E)s{z) + o{p{z)\(E))) (1 + (p(z)ee- 5 °/ £ )) , 



where 



Acr n T v ^ 



cos 



ME) 



Th(E) 

In view of (4.17) and (4.18), we have \(E) = A 7r (l + o(l)). This and (4.35) imply (4.31) and (4.32). 
The first estimate in (4.33) is proved in the same way as the second estimate in (4.25). 
Prove the second estimate in (4.33). As when proving the third estimate in (4.25), one checks that, 
for j 6 {1,2}, \rj\ < Cpe~ 5l l £ and \rj\ < C. Recall that IQ22I < C. These observations and (4.7) 
imply that \R\ < C\p u (z,E) - 1| + \r 2 \ < Cpe~ Sl / £ . 

The "rough" estimates (4.34) follow from the already obtained and (4.26). This completes the proof 
of (4.31) - (4.34). 

Now, assume that e is so small that \g(z,E)\ < 1/2 for all z and E in the case of Proposition 4.1. 
This is possible in view of (4.34). Then, substituting representation (4.31) into (4.28), and taking into 
account (4.33), we get 

v = F(E) + X n s(z) + o(p\ n ) + — — - — — (l + g(z,E)) 

l + g(z,E) 

= [F(E) + X^z) + o{pK) + 0(R(z, E)) + 0(g{z, E))) (1 + g(z, E)) 
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with 

F(E) = F(E) + (Qn(E) + Q 22 (E)). 
In view of (4.32) and (4.33), this implies (4.8). 

Now, we only have to check (4.9) to complete the proof of Proposition 4.1. For sufficiently small e, 
the representation for F in (4.9) follows from 



))■ 



(4.36) F(E) = 4a 7T (T h (E 7T )r 1 cos (^j^) (E - E n ) (1 + o(l 

(4.37) Q n (E) + Q 22 (E) = -2a n A n (V) sin C^j 2 ^) + o(l). 

The formula (4.36) follows from (4.17), (4.14) and (4.18). To prove (4.37), we note that, by (3.11), 

Qu(E) + Q 22 (E) = (9 + 1/9) cos - $ )/e) • 

This in conjunction with (4.15), (4.13), (2.21), (6.3) and (4.3) yields (4.37). 
Finally, the asymptotics for F' in (4.9) follows from 

(4.38) F'(E)=FUl + o(l)), \F^\>Ce- 2 e s ^, \Q' U (E)\ + \Q' 22 {E)\ < Ce~ x . 

Prove the first of these estimates. It follows from Lemma 2.1 and estimates (4.19), (4.13) and (4.16) 
that 

±h\E) e 

Now, using (4.17), (4.15) for v = vr, (4.18) and the estimate & n (E) = & n (E n )(l + o(l)) (following 
from Lemma 2.1), we get 

This and the definition of F n imply the representation for F' in (4.38). The estimate for \F^\ follows 
from the definition of F n and the estimates (4.16), (2.29) and (2.25). The last estimate in (4.38) 
follows from (2.24), (2.32) and the Cauchy estimates for E 9(E). 

This completes the proof of Proposition 4.1. □ 

4.2. The location of the spectrum. We now prove Theorem 1.2. Therefore, we apply Proposi- 
tion 3.1 to the scalar equation with the coefficients p u and v u computed in section 4.1. 
Let Jl the subinterval of J described by (4.2). One has 

Lemma 4.3. The spectrum of H ZtE in Jf is contained in the interval described by 

(4.39) |i^(£)| < (2 + |A*|) (l + o(l)), 
where o(l) is independent of E and E n (satisfying (4.1)j. 

Proof. First, we find r, a subset of J^, where M u , v u and p u satisfy the assumptions of Proposition 3.1. 
Hence, r is in the resolvent set of (0.1). 

Recall that (z, E) i— > p u (z, E) and (z, E) i— > v (z, E) are real analytic as the matrix (z, E) i— > M u (z, E) 
is. Therefore, the equalities ind p u (-,E) = and mdv u (■, E) = automatically follow from the 

inequalities min \p u (z, E)\ > and max \p u (z, E)\ < — ( min \v u (z, E)\ 
zm z m 4 \zgr 

Furthermore, by (4.10), the first of these inequalities is satisfied for all E G J^. So, in J^, the 

assumptions of Proposition 3.1 are satisfied if and only if max \p u (z, E)\ l l 2 < — min \v u (z, E)\. 

zeR 2 zgk 

Corollary 4.1 yields 

(4.40) max \p u (z)\ 1/2 < 1 + o(l), 

(4.41) \ mm \v u {z)\ > (mm \F n (E) + K sin(x)| + o(l)(l + 1^1) 

Z z£K Z yx6K 

where o(l) is independent of E and E w . So, v u and p u satisfy the assumptions of Proposition 3.1 if 
E satisfies the inequality of the form \F n (E)\ > (2 + (A^l) (1 + o(l)), where o(l) is independent of E 
and E n . Now, Proposition 3.1 implies the statement of Lemma 4.3. □ 
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Lemma 4.3 and the definitions of X n and F n , namely (4.4) and (4.5), imply that, in J^, the spectrum 
of H z £ is contained in /„■, the interval described by 



$;(E 7r ) (E - E n ) - A n (V) T h (E n ) tan 



< 2 



T h (E n ) 



cos 



+ T V>V (E V ) (l+o(l)), 



where o(l) depends only on e. The interval /„- is centered at the point 



(4.42) 

and it has the length 
(4.43) 



E w — E v + 



2e 



eA n (V)T h (E n ) 
Th(E n ) 



tan 



K(E.) I 2 



cos 



+ r v , 7r (E 7r ) (i + o(i)). 



This completes the proof of Theorem 1.2. 
Note that 

(4.44) \E n - E w \ + < Cee~ 5o/£ . 

These estimates follow from (4.42), (4.43) and estimates (2.25), (2.29) and (4.16). 
Finally, we note that, using (4.42), one can rewrite (4.4) as 

'$o(En)\ K(E«) 



□ 



(4.45) 



F^E) 



4oy 
Th{E n ] 



cos 



(E — E n ). 



4.3. Computation of the integrated density of states. We now compute the increment of the 
integrated density of states on the intervals described in Theorem 1.2 and, thus, prove Theorem 1.3. 
We use the approach developed in [12]. One has 

Proposition 4.2. Pick two points a < b of the real axis. Let 7 be a continuous curve in C+ connecting 
a and b. 

Assume that, for all E € 7, one can construct a consistent basis such that the corresponding mon- 
odromy matrix is continuous in (z,E) £ M. x 7 and satisfies the conditions 

2 



(4.46) 



min|Mi2| > 0, max|p(z)| < I — min \v{z 



z til 



1 



where p and v are defined by (3.1) with h from (2.5). Assume in addition that the coefficients of M 
are real for real E and z. Then, one has 



(4.47) 



N(b) - N(a) = - ^ j ar S v ( z > E) dz 



where /| 7 denotes the increment of f when going from a to b along 7. 

Proof. In [12], we proved a more general result; we assumed that, for all (z, E) G 8x7, the monodromy 
matrix satisfies the conditions of Lemma 3.1 and got the formula 



(4.48) 

where G is the continued fraction 
(4.49) G{z)=v{z) 



N(b) - N(a) = -^J arg G(z, E) 



dz 



v (z — h) 



p(z ~ h) 



p(z - 2h) 



v(z- 2h) 

Such continued fractions were studied in [4]. It was proved that, if the functions z \—* p(z) and z *— > v{z) 
are continuous and 1-periodic and if they satisfy the conditions (3.3), then 

• the continued fraction z 1— > G{z) converges to a continuous 1-periodic function uniformly in M; 

• if p and v depend on a parameter E, if they are continuous in (z, E) in some domain D, and 
if, for all (z, E) S D, they satisfy conditions (3.3), then (z, E) 1— ► G(z, E) is also continuous in 
D. 
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i in 



• for z G M, one has 

(4.50) \G(z) - v{z)\ < £ min \v(z)\ - \l I £ min ) - max |p(z) 

z 2gR y yz ^gk y zgM 

Now, turn to the proof of (4.47). As, in our case, v(z,E) and p(z,E) are real for real z and E, we 
conclude that (1) indv = indp = (which follows from (4.46)); (2) the right hand sides in both (4.47) 
and (4.48) belong to e/2-7rZ. The first observation and (4.46) imply that, for all (z,E) G K x 7, 
the monodromy matrix satisfies the conditions of Lemma 3.1. In view of the second observation, 
formula (4.47) follows from (4.48), the continuity of (z,E) 1— > G(z, E)/v u (z, E) and the inequality 

sup 2gK ^ Z, ^u ^ E)\ E ^ ^ 1 van d f° r an E G 7. And, the last one follows from (4.50) and the second 
condition from (4.46): 

\G(z,E)-v u (z,E)\ < 2 m a x zm \p u (z,E)\ < 1 
2GRj e ^ \v u (z,E)\ mm zm \v u (z,E)\ 2 

This completes the proof of Proposition 4.2. □ 

4.3.1. The computation. Let E n be as in the beginning of section 4 and, in particular, be such that (4.1) 
is satisfied. As above, let J%, be the subinterval of J described by (4.2). 

As seen in the previous section, in JJ, the spectrum of H ze is contained in I n , the interval centered 
at E w (see (4.42)) of length \I„\ (see (4.43)). 

To compute the increment of the integrated density of states on I n , we use Proposition 4.2 and choose: 



7 



$E G C+ : \E- E n \ = l e a e -*/e| 



Let a < b be the ends of 7. Then, by (4.44), one has C (a, b). We prove 

Lemma 4.4. On 7, the monodromy matrix M u and the functions p u and v u satisfy the condi- 
tions (4.46). 

Recall that the integrated density of states of H Z £ is constant outside the spectrum of H Z £ . So, its 
increment on I n is equal to its increment between the ends of the semi-circle 7. And, in view of 
Lemma 4.4, the latter is given by the formula (4.47). In view of this formula, to prove Theorem 1.3, 

it suffices to check that j} &rgv u (z, E) dz = — n. This follows from 

7 

Lemma 4.5. For (z,E) G M x 7, one has 

(4.51) v u (z,E)=F w (E)(l + o(l)). 

Indeed, note that for z G M and E G M, the functions F n and v u take real values. Therefore, the 
estimate of Lemma 4.5 implies that J &rgv u (z, E) dz = a,rg F n (E)\ ry . In view of (4.45), the last 

quantity is equal to (E — i^r)| = — vr. So, to complete the proof of Theorem 1.3, we have only to 



check Lemmas 4.4 and 4.5. They will follow from 

Lemma 4.6. For (z, E) G R x 7, one has 

(4.52) 1^(^)1 > Ce a ~ 2 . 

Proof. The lower bound for ji^i?)! follows from (4.45), the definition of 7 and the estimates (2.29), 
(4.16) and (2.25). □ 



Proof of Lemmas 4-5. Prove the asymptotic representation for v 11 . Therefore, we first derive an upper 
bound for the ratio K/F^E). By (4.5) and (4.45), we get |A 7r /F 7r (E)| = ,,f;^\ , . Now, the 
definition of 7 and the estimates (2.29) and (2.25) imply that 

(4.53) \K/F^{E)\<Ce 1 - a e- So ' £ . 

So, the ratio is small when e tends to 0. The representation (4.51) follows from (4.11), (4.53) and (4.52). 
This completes the proof of Lemma 4.6. □ 
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Proof of Lemmas 4-4- By Proposition 4.1, for sufficiently small e, for Be 7 and z G R, one has 
M^ 2 (z, E) / 0. Finally, for sufficiently small e, for all £ e 7, from (4.52) and (4.7), it follows that 

max 1 0(2, < — min E)\ 2 . This completes the proof of Lemma 4.4. □ 



4.4. Computation of the Lyapunov exponent. We now derive the asymptotics of the Lyapunov 
on the interval I n , i.e., prove formula (1.19), and, thus, prove Theorem 1.4. 



4.4.1. Preliminaries. To compute Q(E,e), we use Theorem 2.1 and compute the Lyapunov exponent 
of the matrix cocycle defined by the monodromy matrix M U (-,E). It appears to be difficult to compute 
directly 9(M U '(•, E),h): one can obtain only rough results. However, using the scalar equation with 
the coefficients v u and p , one can construct another matrix cocycle that has the same Lyapunov 
exponent as (M u (-,E),h) and for which the computations become much simpler. 



4.4.2. The Lyapunov exponent and the scalar equation. In this section, we assume z 1— ► M(z) to be a 
1-periodic, SL(2, R)-valued, bounded measurable function of the real variable z. Let h is a positive 
irrational number. We check the following simple 

Lemma 4.7. Assume that there exists A > 1 such that 

(4.54) Vz G R, A" 1 < Afi 2 (s) < A. 

In terms of M and h, construct v and p by formulae (3.1). Set 



(4.55) 



N / z \ = ( v ( z )/y/f>(z) -y/rf/i 
1 ; V Vy/M 



Then, the Lyapunov exponents for the matrix cocycles (Af , h) and (N, h) are related by the formula 
(4.56) 9{M,h) = 6(N,h). 



Proof. Let 



One has 



H(z) = ^_(M 12 (z) 
H(Z) M l2 {z) \M 22 (z) -1 



(4.57) M(z) = e l{z) H(z)N{z)H~ l (z - h), l(z) = - logp(z). 

Note that, under the condition (4.54), 

\l(z)\ < log A < 00, Vz G R, 
and that l(z) is 1-periodic. As h is irrational, by Birkhoff's Ergodic Theorem ([23]), one has 



(4.58) hm \y^l( z + jh)= f 



l{z)dz 



for almost all z G R. As 2l(z) = log p(z) = logMi 2 (z) — logMi 2 (z — h), the integral in (4.58) 
vanishes. This, the definition of the Lyapunov exponent (2.7), relation (4.57) imply relation (4.56). 
This completes the proof of Lemma 4.7. □ 

Now, for p = p u and v = v u , we construct N u by formula (4.55). Relations (2.8) and (4.56) imply 
that the Lyapunov exponent for the operator H z e is given by the formula 

(4.59) Q(E,e) = ^6(N u (;E),h). 

In the next two subsections, we prove a lower and an upper bound for 9(N U (■, E),h). They will 
coincide up to error terms, and, thus, yield the asymptotic formula for Q(E,e). 
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4.4.3. The lower bound for the Lyapunov exponent. Here, we prove that, in the case of Theorem 1.2, 
for E E I n , the Lyapunov exponent admits the lower bound: 

(4.60) 6(N u (-,E),h) > log + \X n \ + o(l). 
Therefore, we use the following construction. 

Assume that a matrix function M : C — > SL(2, C) is 1-periodic in z and depends on a parameter 
e > 0. One has 

Proposition 4.3. Let eq > 0. Assume that there exist yo and y\ satisfying the inequalities < yo < 
yi < oo and such that, for any e E (0,£o) one has 

• the function z — ► M(z,e) is analytic in the strip {z E C; < Imz < yi/e}; 

• in i/ie sirip {z E C; yo/ e < Imz < yi/e}, M(z,e) admits the following uniform in z represen- 
tation 

(4.61) M(z,s) = X(e)e 2mmz -((l °_) + o(l)) , e - 0, 



where A(e) and m are constant; m is integer (independent of e). 
Then, there exit a E\ > smc/i i/iai, i/ < e < £\ , one has 

(4.62) 0(M,fr)>log|A(e)|+o(l); 

the number £\ and the error estimate in (4.62) depend only on eq, yo, y\ and the norm of the term 
o(l) in (4.61). 

This proposition immediately follows from Proposition 10.1 from [12]. Note that the proof of the latter 
is based on the ideas of [25] generalizing Herman's argument [16]. 

We apply Proposition 4.3 to the matrix N u (z, E). Therefore, we fix 2/2 and y\ so that < y2 < y\ < 
y < 5i/(2ir), where 5\ is the constant from the Proposition 4.1. Then, the estimate (4.60) follows from 
Proposition 4.3 and 

Lemma 4.8. Assume that X w > 1. In the strip y^ < Imz < y\, for E E 1^, the functions p u 
satisfies (4.10) and v u admit the asymptotics: 

(4.63) v u (z, E) = K e- 2 *^-**^)) (1 + G (i)). 
These asymptotics are uniform in E n (satisfying (4A)), E and z. 

We postpone the proof of this lemma and complete the proof of the estimate (4.60). If jA^I < 1, the 
estimate (4.60) gives a trivial lower bound as the Lyapunov exponent is always non-negative. So, it 
suffices to prove (4.60) in the case | | > 1. Substituting (4.10) and (4.63) into (4.55), for E E I n and 
yijz < hn z < yi/e, one obtains 

" A N 



0<+°« 



N u {z) = X n e -2Ti(z-^(^)) 

as z n is real and \e 2mz \ > e 27ryi ^ £ > 1. Proposition 4.3 then implies (4.60). □ 

Proof of Lemma The first statement is taken from Corollary 4.1. Let us prove (4.63). First, we 
recall that, as E E I n , one has (4.39). On the other hand, for Imz > 2/2/^ > 0, one has 

(4.64) X w sin(27r(z — h — z 7r (E 7T ))) = ^ e -2*i(*-h-*0 (1 + o(1)) . 

Note that, as | A^l > 1 and z n E M, the right hand side is exponentially large as e — > 0. Then, in the 
strip 2/2 /£ < Imz < yi/e, for E E (4.11), (4.39) and (4.64) imply (4.63). This completes the proof 
of Lemma 4.8. □ 

4.4.4. The upper bound for the Lyapunov exponent. We now prove that, in the case of Theorem 1.2, 
the Lyapunov exponent admits the upper bound 

(4.65) 6{N u (-,E),h) < log + |A^| +C. 

This upper bound follows from the definition of Lyapunov exponent for matrix cocycles (2.7) and the 
estimate 

sup sup IliV^z,^)!! < C{\X n \ + l), 
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which follows from (4.55), (4.10) and the estimate 



sup sup \v u (z, E)\ < C(\X n \ + 1), 

which follows from (4.11) and (4.39). This completes the proof of (4.65). □ 

4.4.5. Completing the proof of Theorem Estimates (4.60) and (4.65) together with the represen- 
tation (4.59) imply the uniform representation 

V£e 4, Q(E,e) = ^log+|A 7r (E 7r )| + 0(e). 

In view of (4.5), to complete the proof of Theorem 1.4, it suffices to check that 



cos 



- inf|^-4 |, 

£ I 



which follows from the definition of the points (Eq^)i and from (2.24) and (2.25). □ 

4.5. Absolutely continuous spectrum. We now turn to the proof of Theorem 1.5. The idea is to 

find a subset of I n where E i— > @(E,e) vanishes. Then, by the Ishii-Pastur-Kotani Theorem ([23]), 

this subset is contained in the absolutely continuous spectrum of the ergodic family (0.1). 

As before, we assume that h is defined by (2.5), and that the functions p u and v u are the coefficients 

of the scalar equation equivalent to the monodromy equation with the matrix M u . 

As in the previous subsection, to analyze Q(E, e), we use the matrix cocycle (N u (-, E),h), the matrix 

N u being defined by (4.55) for M = M u . Recall that @(E,e) is related to 6(N U '(•, E), h), the 

Lyapunov exponent of this cocycle, by the formula (4.59). 

First, under the conditions of Theorem 1.5, we check that, up to error terms, N u is independent of z. 
This allows then to characterize the subset of where 6{N U , h) = by means of a standard KAM 
construction found in [12]. 

4.5.1. The asymptotic behavior of the matrix N u . We need to control the behavior of the matrix N u 
for bounded |Imz| and E near the interval /„-. One has 

Lemma 4.9. Fix c > 0, x > and r > 0. For e sufficiently small, if E n satisfies (4.1), and if 
e log X n (E 7T ) < — c, then 

N u (z,E) = N (E)+N 1 (z,E), 

(4.66) No{E) = fHE) -l\ gup \\ Nl ( Zi E)\\ <Ce-^, 

\ J \E-E^\<x%\ 

|Imz|<r 

where the constant n is defined by rj = min{c, 5±}, and F is the function from (4.8). 

Proof. It suffices to prove, that under the conditions of the lemma, there exists C > such that, for 
e sufficiently small, one has 

(4.67) sup \p U (z,E) - 1| < Ce-^ £ , sup \v u (z, E) - F(E)\ < Ce^ £ , 
\E-E 7I \<x\I 7I \ \E-En\<*]I*\ 

\Imz\<r \lmz\<r 

(4.68) sup \F(E)\ < C. 

\E-K\<x%\ 

Begin with the proof of (4.68). Recall that, for E in the ghborhood of E n , one has (4.9). 

On the other hand, the interval I n is located in the (Cee _<5 °/ £ )-neighborhood of E^, see (4.44). So, it 
suffices to prove (4.68) with F replaced by F w . 

Recall that I n is centered at E n , see (4.42), and that, by (4.45), one has F^^^) = 0. The esti- 
mate (4.39) is an estimate for F V (E) on the interval As E i— > F n (E) is affine, it implies that 

sup \F n (E)\ <x(l + |A ff |)(l + o(l)). 

\E-E 7I \<x\I 7T \ 
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As \X W \ < 1, this implies (4.68). 

Let us prove (4.67). The representation (4.8) and estimate (4.68) imply that, for some C > 0, 

sup \v u (z,E) - F(E)\ < Cee- 5o/£ sup \F(E)\ + C K + C e~ 5l/£ < 
\E-E v \<>t\I*\ \E-E*\<x\I v \ 
|Imz|<r 

< C{ee- 5 » /e + e~ c/e + e~ Sl/£ ). 

In view of (4.6) and the definition of r], this expression is bounded by C e~ v ^ £ . This proves the second 
estimate from (4.67). The first one follows from (4.7), (4.6) and the definition of rj. Lemma 4.9 is 
proved. □ 

4.5.2. The KAM theory construction. Here, we formulate a corollary from the construction developed 
in section 11 of [12] that is based on standard ideas of KAM theory (see [5, 2]). 

Let / C K be a bounded interval. Fix r > 0. Let S r be the strip {z £ C; |Imz| < r}. We consider A, 
the set of 2 x 2-matrix valued functions (z, <p) G S r x 1 1— > M(z, p) that are 

(1) analytic and 1-periodic in z £ S r ; 

(2) analytic in p in V(I), a complex neighborhood of /; 
a b 
b* a*. 



(3) of the form 



(e lip \ 

Let D = I _j , and let A £ A satisfy X(A) = sup \\A(z, <p)\\ < oo. 

V u e J <pev(i),zes r 

Fix < h < 1. For z61h ^{z) G C 2 , a vector function, consider the equation 

(4.69) ip{z + h) = (D + A)(z)ip(z). 
Define 

H(u) := {h e (0, 1); min \h - l/k\ > n/k 3 for k = 1, 2, 3 ... }. 

One has 

Proposition 4.4. Fix a S (0,1). There exists Ao(r, a, I) > smc/i £/ia£, /or any ^4, D and h chosen 
as above and satisfying the conditions 

(1) det(D + A) = 1, 

(2) A = X(A) < X (r,a,I), 

(3) h e H(X a ) 

there exists C I, a Borel set of Lebesgue measure smaller than A CT / 2 and such that, for all tp G 
I\ $oo? equation (4.69) has two linearly independent bounded solutions. 

This proposition immediately follows from Proposition 11.1 of [12]. The constant Xo(r,o~,I) depends 
only on the length of /, but not of its position. 

Proposition 4.4 implies 

Corollary 4.2. In the case of Proposition 4-4> f or all p £ I \ &oo, the Lyapunov exponent of the 
cocycle (D + A, h) is zero. 

Proof. Let ^(z) be the matrix the columns of which are the vector solutions defined in Proposition 4.4. 
Then, ^>(z) is a matrix solution of (4.69). As the vector solutions are linearly independent, det ^?(z) ^ 
for all z£R. For I £ Z, put x(0 = + Ih). Then, X (Z + 1) = (D + A)(hl + z)xQ), and, as #(z) is 
bounded, for L > 1, we have 

|| (D + A)(Lh + z) ■ ■ ■ (D + A)(h + z)(D + A)(z)\\ = \\ X {L + l)x -1 (0)|| < C. 

Now, the statement of the corollary follows from (2.7), the definition of the Lyapunov exponent. □ 

4.5.3. The proof of Theorem 1.5. The idea is the following. Let S be a constant matrix such that 
det 5/0. Clearly, 

(4.70) 0{N U , h) = e(S^N u S, h). 

Recall that N u admits the representation (4.66). We shall choose S so that the matrices 

(4.71) D = S- 1 N S and A = S- 1 N 1 S 
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satisfy the assumptions Proposition 4.4. Then, we apply Corollary 4.2 to the so constructed matrix 
D + A. We divide the analysis into "elementary" steps. 

Diagonalization. Let F° be a point of I n such that 

-1 < F(F) < 1. 

Then, in V°, a neighborhood of F°, one can define an analytic branch of the function E \— ► y>(F) 
solution to 

(4.72) cos ip(E) = F{E). 

In V°, the exponentials e ±J</3 ^^ are the eigenvalues of the matrix Nq(E) (see (4.66)); the columns of 
the matrix 



S(E) 



e i<p(E) e -ME) 



1 1 

are its eigenvectors. Define D and A by (4.71). Clearly, 

(4-73) D W=[o e X 

As E i— > N\{E) is real analytic, -A(z, F) has the form 

A 

For some C > 0, one has 



a 6 

6* a* 



2|I m <p(E)| 

(4.74) VFeF , su V \\A(z,E)\\ <C— -— sup ||JVi(z, E)\\. 

z m \sm ip(E)\ zeR 

A change of variables: E — > <p. Now, we change the variable E to <p, and check that, as a function of 
99, A satisfies the conditions of Proposition 4.4 and Corollary 4.2. We use 

Lemma 4.10. Fix x < 1. There exists £q > such that, for < e < £0 the following holds. Let E n 
satisfy (4.1). Let I C R be the interval centered at E^ and of length x\I n \. Then, 

• in a neighborhood of I, there exists a real analytic branch 0/ 99(F); it is monotonous on I; 

• there exists a positive A = A(x) such that <p(I) C (A,7r — A); 

• (p 1 — ► E{ip), the function inverse to E 1— > v?(F) is analytic in V{L), the A/2-neighborhood of 
the interval ip{I), and maps V(I) into the (C\I n \) -neighborhood of I n . 

As F(F) is real analytic, Lemma 4.10 immediately follows from (4.72) and 

Lemma 4.11. Fix >c\ G (0, 1). For e sufficiently small, the following holds. Let E n satisfy (4.1) and 

define B = {E G C; |F - F^l < ± xri 1 4 1 } . 

Then, F bijectively maps B onto F{B), and one has 

(4.75) sup \F(E)\ < Xi + o(l), and, /or |F - F^| = ^ IFJ, |F(F)| = x x + o(l). 

Proof. Fix < a < 1. By (4.44), F is contained in the e a e _5o//£ -neighborhood of E w . Therefore, 
F'(F) admits the representation (4.9). This implies that F is a bijection of B onto F{B). Indeed, 
assume that, in B there exist E\ and E2 such that E\ 7^ F2 and F(E±) = F^E?). Then, one has 

= F(F 2 ) - F(F X ) = [ 2 F'(E)dE = F' n [ \l + o(l))dE = F^(F 2 - Fi)(l + o(l)) ^ 0. 

So, we get a contradiction, and F is a bijection. 
Estimates (4.75) follow from the following facts: 

(1) the representation for F from (4.9) holds on B (as B is contained in the ghborhood 
ofFV); 

(2) E 1 — ► F n (E) is affine, and vanishes at F^, the center of /„•; 

(3) at the ends of one has \F n (E)\ = 1 + o(l) (by (4.39), which is the definition of I n , and as 
K = 0(e-^)). 
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This completes the proof of Lemma 4.11. □ 

Now, turn to the matrices D and A defined by (4.71). Make the change of variables E — > <p so that 
E = E(ip). Consider these matrices as functions of <p in V(I). Then, for e sufficiently small, they 
satisfy the conditions of section 4.5.2: 

• z i — y A(z,ip) is analytic and 1-periodic in S r (as z i— > N u (z,E) is analytic in the strip 
{|Imz| < y}); 

• (p i— ► A(z, ip) is analytic in V(T) (as p i— ► E(<p) is analytic in V(I), ip(V(I)) is in the (C| | )- 
neighborhood of E^, and as i? i— > N u (z, E) is analytic in this neighborhood); 

• A has the form ^ j (as <p \— > -E(y) is real analytic, and as E 1 ^4(z, £?) already had this 

form) ; 

• D is given by (4.73); 

• \(A) < ^e-v/z (by (4.66), (4.74) and Lemma 4.10); 

• det(D + A) = 1 as D + A = 5 _1 iV t/ 5 and detiV^ = 1 by (4.55). 

The Diophantine condition on e. To apply Corollary 4.2, we have to impose a Diophantine condition 
on the number 2tt/e. Fix two positive numbers a and b. Consider the set 

D(a,b) := lee (0,1) : min — -l/k > ^e~ b/e , k = 1,2,3... 1. 

{ leN s k 6 J 

It can be easily checked 

. . mes (D(a, b) n (0, e)) ( h u\ 

(4.76) v v y — \-LJL = i + o [e 1 ) when e -> 0. 

The derivation of (4.76) is similar to the estimates in section 4.4.6 of [12]. 

Fix < a < 1. For e £ D((C/ A) a , err/), the number /i defined by (2.5) belongs to the class H(/i) with 

Completing the proof of Theorem 1.5. Let A and D be as constructed above and £ £ D. Then, for 
the matrix cocycle (D + A, h), the conditions of Corollary 4.9 are satisfied provided e is sufficiently 
small. So, for e sufficiently small, there exists $oo, a subset of I of measure uniformly small with 
A(A) < £ e _r?//£ , such that, for all 99 € I \ $00, the Lyapunov exponent 0{D + A, h) is zero. 
By (4.70) and (4.59), this implies that @(E), the Lyapunov exponent for the family of equations (2.1), 

f dE 

is zero on ip(I) C 1^ outside a set of Lebesgue measure m := / ——dip. 

7 $00 dip 



dE 



dtp 



(<p) < C\I n \ for ip £ <p(I). So, m = o(\I\) 



The Cauchy estimates and Lemma 4.10 imply that 
where |/| denotes the length of /. 

As m is small with respect to |/| and as h in the definition of / can be chosen arbitrarily close to 1, 
we conclude that @(E, e) is zero on 1^ outside a set the measure of which becomes small with respect 
to l-Z^I as e tends to zero in D(rj). This completes the proof of Theorem 1.5. □ 



5. Computing the monodromy matrices 

In this section, we prove Theorem 2.2. As we have seen, to study the spectrum of (0.1), one has 
to compute the coefficients of the monodromy matrix up to terms that are exponentially small (in e) 
whereas these coefficients are exponentially large outside small "resonant" neighborhoods (where the 
points {E w (l)}i are exponentially close to {Eq(1')}i/). To achieve such an accuracy, we use a natural 
factorization of the monodromy matrix into the product of two simple "transition" matrices and carry 
out a rather delicate analysis of the properties of their coefficients. 

Below, we always work in terms of the variables 

(5.1) x := x — z, and £ = ez. 
In these variables, equation (2.1) takes the form 

d 2 

(5.2) — r^V'(^) + (V(x) +acos(ex + C))i(j(x) = EMx), x£R, 

dx A 

30 



The advantage of the new variables is that now we can study solutions of (5.2) analytic in (. 
In terms of variables (5.1), the consistency condition (2.2) takes the form 



^i0, C) 



(5.3) ij> j (x + l,Q = il> j (x,( + e). 
The definition of the monodromy matrix, (2.3), turns into 

(5.4) V(x,( + 2ir)=M((,E)y (x,C), *(x,() 

and, now, the monodromy matrix is e-periodic: 

M(C + e,E) = M(C,E), VC. 

5.1. Transition matrices. Here, we describe the factorization and the asymptotics of the transition 
matrices. 

5.1.1. Factorization. Here, we describe a natural factorization of the monodromy matrix under the 
assumption (TIBM). 

Two consistent bases. In section 7, we pick a point E* in J and show the existence of 14, a neigh- 
borhood of E*, such that, for E G V*, there exists two consistent bases which will be indexed by v in 
{0, 7r}. Let us describe some properties of these bases; they will be central objects in this section. 
Fix v G {0, 7r}. The corresponding basis consists of two consistent solutions to (5.2), say (x,(,E) 
f u (x,(,E) and (x,(iE) h-> f*(x, (, E); the second solution is related to the first one by the transfor- 
mation (2.10). For any the function (£, E) i— ► f v (x, (, E) is analytic in the domain 

(5.5) {( £ C : |ImC| <Y}x V*, 

where Y satisfies the inequality Y > Ym (recall that Ym is defined in (2.13)). 

Definitions of the transition matrices. As both pairs ({fu, /^}W{o,tt} are bases of the space of 
solutions of (5.2), one can write 



(5.6) F n (x,C + 27r,E)=T n (C,E)F (x,C,E), F (x, £, E) = T ((, E) F n (x, (, E), F u 



For v G {0,7r}, the 2 x 2-matrix valued function (C,E) i— > T U ((^,E) is independent of x. We call it a 
transition matrix. 

Discuss the basic properties of a transition matrix. As the basis {f u , /*} is consistent, for all E, 
£ i — ^ T U (£,E) is e-periodic. It is analytic in the domain (5.5). Finally, as the consistent solutions 
f v and /* are related by the transformation (2.10), T v enjoys the same symmetry property as the 
monodromy matrix (see (2.9)); we write 



a v b v 
b* a* 



Factorization of the monodromy matrices. For v G {0,7r}, we denote by M v the monodromy matrix 
corresponding to the base {f u , /*}. The definitions (5.4) and (5.6) imply that 

(5.7) M n (C) = T n (() T (C), M (C) = T (C + 2vr) T^C). 

Clearly, the monodromy matrices share the basic properties of the transition matrices: they are e- 
periodic in £, analytic in the domain (5.5) and have the form (2.9). 

Note that, once transformed back to the z-variable, the monodromy matrices are analytic in the 
domain {C G C; |ImC| < Y/e} x V*. 
Finally, by (2.4) and (5.7), one has 

(5.8) detT detT v = 1. 

The motivation for considering the factorizations is the following. The solutions /o and f n are con- 
structed so that fo has a simple asymptotic behavior in the strip {— ir < Re(" < tt}, and f n has a 
simple asymptotic behavior in the strip {0 < Re£ < 2tt}. In result, formulae (5.7) give factorizations 
of the monodromy matrices in terms of factors with simple asymptotic behavior. 
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5.1.2. Asymptotics of the transition matrices. We now describe the asymptotics of the transition 
matrices (T v ) ve i Q ^y. Therefore, we shall use the conventions introduced in (2.11), (2.12) and (2.13) 
in section 2.2. We need a few more notations. 

1. Asymptotic notations. We shall use all the notations introduced in section 1.10. 

2. "Analytic" notations. Pick zq G K and let Vjj be a complex neighborhood of zq. Let z i— ► a(z) 
be an analytic function defined and non vanishing in Vq. In Vq, we define two real analytic functions 
z i — > I a I (z) and z i— > 93(a) (z) by 

a(z) = I a I (2) exp(i<£>(a)(z)) 

such that I a I (z) = \a(z)\, and ip{a)(z) = arga(z) when z G Vo n E. 

5. "Fourier expansion" notations. The transition matrices being e-periodic, we represent their Fourier 
expansion in the form 

(5.9) a„(C) = a„,-i(C) + a,, + a„,i e 2 </ £ + a„ )2 (C), MO = 6„,-i(C) + k,o + e 2 " C/£ + MO, 

where we single out the sum of Fourier terms with negative index, the zeroth and the first Fourier 
terms and the sums of Fourier series terms with index greater than 1. 

One has 

Theorem 5.1. Pick E* G J. There exists V* ; a complex neighborhood of E*, and Y > Ym such that, 
for sufficiently small e and v G {0,7r}, there exists {fi>,f„}, a consistent basis of solutions to (5.2), 
having the following properties: 

• the basis {f v ,f„} and the transition matrices T v are defined and analytic in the domain (5.5); 

• the determinant of T v is independent of £ and e; it is a non-vanishing analytic function of 

E g K; 



• one has 



(5.10) 



ant 



(5.11) 



(5.12) 



I a ufi I = exp \^Sh, v + 0(1) \ , I &„ j0 \ = exp + 0(l)j , 

I a„ 7 i I = exp (^(Sh,u ~ S VjV ) + O(l) J , I b Vj i \ = exp \^(Sh,v - S v>u ) + 0(1)J , 

V>M = ~y £ ( $ o - + 0(1), v(6o,i) = ~Y e ($0 + $7r) + 0(1) ' 

¥>K,l) = (^ - ^0 - 4vr 2 ) + O(l), y?(^,i) = + $ - 4vr 2 ) + 0(1), 

where 0(1) denotes functions real onV*C\M. and analytic in E G V*; 
• moreover, 

(5.13) a„_i(C) = o(a„ i0 ), V-i(C) = 0(6^0), 0^2(0 = o(p{(/e) a V)1 ), b Vj2 (() = °(p((/ £ ) Vi)- 
All the above estimates are uniform in E and £ in the domain (5.5). 

Theorem 5.1 is proved in sections 7—11. 

When studying the spectral properties of (0.1), we always assume that E satisfies 

(5.14) E G := V* Pi {|ImE| < e}. 
One proves 

Corollary 5.1. Pick u G {0,7r}. For sufficiently small e, in the case of Theorem 5.1, for E G V+, 
one has 

/r ir\ I I 1 It I 1 I I \tl',v\ 1, I \tv,v\ 

I 5 - 15 ) K,0| - 77 1> \0u,0\ - 77 f, 7: — f) 1^,1 1 



\th,uY \th,u\' \th,vY \th,i 



■V 

where all the tunneling coefficients are computed at the point Re-E instead of E. 
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Proof. The functions E \— > Sh,o(E), E \— > Sh,n(E), E ^ &o(E) and E \— > & n (E) are independent of e 
and analytic in a neighborhood of J. So, for sufficiently small e, for E E V^f, one has 



= j*,/(Re£;)/e| 



IW(^)| - \td,u(ReE)\, \e~ 

for E {0, 7r} and for (i E As the phase integrals are real analytic, one has | e i ^^(^)/ £ | x i. 

Estimates (5.15) follow from these observations and representations (5.10) — (5.12). This completes 
the proof of Corollary 5.1. □ 

5.2. Relations between the coefficients a v and b u of the matrix T v . It appears that, with 
a great accuracy, the coefficients a u and b v are proportional. This makes the factorizations (5.7) 
extremely effective. Recall that Y m is defined in (2.13). Define 

b„((,E) 



(5.16) 
One has 



RACE) 



a„({,E) 



Proposition 5.1. Pick v E {0,7r}. Fix < y < Y m . For e sufficiently small, in the case of 
Theorem 5.1, for |Im£| < y and E E V£ one has 



(5.17) 



Rv((,E)=e 



i{ip{b v fi)-ip{a v fl)) 



1 



det T v 



where 

(5.18) O v = 0(t\ v , T YP (C/e), t 2 htV t v , v p{C/e)). 

Proof. In this proof, we assume that E E V*. We set 



and note that 

(5.19) 

and 

(5.20) 



Y VtV (E) = ^-S v , v (ReE), 

Z7T 



< y < Y m < Y VtV (E) < Y M < Y, 



\t v ^{E)\^e-^ Y ^l £ . 



The plan of the proof is the following. We first prove that, for |Im£| < y, 



(5.21) 



Ry 



l + ( e -2,iY-\lr a C\)/e^2,\l m C\/e t ^ t 2^ 



where r v is independent of £ and r„xl. Then, we compute r v with high enough accuracy: we prove 
that 



(5.22) 



1 



det T„ 



2a u0 a 



+ 0(tt e~ 2nY/£ t tl 



;/.0 



Representations (5.21) and (5.22) imply Proposition 5.1. Indeed, to get (5.17), one has to substi- 
tute (5.22) into (5.21) and to take into account that, in (5.22), the second and the third terms in the 
square brackets are bounded by a constant independent of E, £ and e. Note that, from the second 
point of Theorem 5.1 and estimates from Corollary 5.1, it follows that 

(5.23) <Ctl u (ReE). 

To prove (5.21), we use the following observation. 

Lemma 5.1. Pick v E {0,7r}. For sufficiently small e, in the case of Theorem 5.1, one has 

• in the strip |Im£| < Y , each of the functions ( t— > a v (Q,E) and £ i— > b v {C,,E) has one zero per 
period; 

• the imaginary part of the zeros have the asymptotics —Y VyU (E) + 0{e); 

• for any zero of a v , there exists a unique zero of b u such that the distance between them is 
bounded by C et\ u (ReE). 
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We prove this lemma later. In view of the first point of Lemma 5.1, we can represent R v in the form 

27Ti(C-f 6 )/£ _ 1 

(5-24) R„(0=ll„(()p u (0 where n„(C) = e2nl{( _ Ca)/£ _ 1 , 

where ( a (resp. Cf>) 1S one of the zeros of a (resp. b) in the strip {|Im£| < Y}, and is a e-periodic 
function analytic in this strip. The representation (5.21) then follows from the representations: 

(5.25) Ii u {Q = l + 0{e- 2 ^ e t VyU tl v ) for |ImC|<y, 

(5.26) Pu«) = Pv,0 + O(e- 2 ^ Y -\ lm ^ £ ) for |Im£| < Y and p, xl, 

where p^o is the 0-th Fourier coefficient of p. Indeed, to get (5.21), one has just to substitute (5.25) 
and (5.26) into (5.24) and to take into account the fact that the error term in (5.26) is uniformly small 
when |Im£| < y. And the latter follows from (5.19). 

Check (5.25). In view of the second and the third points of Lemma 5.1, and (5.19), for sufficiently 
small e and |Im£| < y, we get 

IMC) - 1 = e^-W/' 1 "f! ( t" Cn)/ . e = O(e-^^y^j) = ofe-^^J) , 



e 27ri(C-Ca)/e _ 1 

where, at the last step, we have used (5.20). This proves (5.25). 

Recall that p v Q be the zeroth Fourier coefficient of p. To prove (5.26), it suffices to check that 

(5.27) |p(C) -p v ,o\ < Ce- 27ry/£ e 2 ^ |ImCI/e for \Im(\<Y and p„, xl. 
Both these estimates follow from the representations 

(5.28) p(() = ^ (l + o(l)) for ImC = -y, p(() = — (1 + o(l)) for ImC = Y. 



Indeed, in view of Corollary 5.1, one has 







a v,0 





>c 1. Therefore, any of the representations (5.28) 
Im£| = Y, we have \p(()\ < C. This bound and 



implies that p v $ x 1; (5.28) also implies that, for 
general properties of periodic analytic functions imply (5.27). So, to complete the proof of (5.26), we 
need only to check (5.28). 

We check only the first of the representations (5.28); the other one is proved similarly. First, we note 
that, for sufficiently small e and ImC = — Y, 

e 27ri(C-Ci,)/e _ 1 

n ^ = e MH a )/e-l =1 + ^ 

Indeed, this follows from the last two points of Lemma 5.1 and (5.19). Now, in view of (5.24), it 
suffices to check that, for ImC = —Y, 

^(0 = ^1(1 + 0(1)), 

which follows from 

aAO = a uA e 2m Hl + o(l)) and b v (Q = b uA e 2m <(l + o(l)). 

We prove only the first one; the second is proved similarly. By Theorem 5.1, Corollary 5.1 and (5.20), 
for ImC = —Y and E € V£, we have 

a v (Q = a^e 2 ^ (l + o(l) + O (^e~ 2 ^\\ = a^e 2 ^ (l + o(l) + O ((VrV 27 ^)) 

= a,,ie 2 ^ (l + o(l) + O {e~ 2 ^ Y ~ Y ^/ e y^ = a v ^ iC > (1 + o(l)) , 

where we have used (5.19). This completes the proof of (5.26) and, thus the proof of (5.21). 
Now, we compute the constant r v from (5.21). First, we prove that 

(5.29) r v rt = 1 - + 0(tV, e^ Yle \ 

This relation follows from the relations 

(5.30) R u R* v = l 



det T v 



a u al 
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and from the fact that, for ( S M, 

(5.31) o„(C) = a„, (l + 0{t VjU )). 

Indeed, recall that all the functions we work with are e-periodic; substituting (5.21) and (5.31) 
into (5.30) and integrating along R over a period, we get 

r„r*(l + 0(e- 2 ^ y / £ , t v ,„ t\ „)) = 1 - _^L(i + (t v , v )). 

a u ,o a v0 

In view of (5.23), this immediately implies (5.29). So, to complete the proof of (5.29), we have only to 
prove the relations (5.30) and (5.31). The relation (5.30) follows from the equalities det T v = a y a* u —b u b* v 
and (5.16). To prove the relation (5.31), we rewrite (5.9) in the form 



(5.32) a u = a u fl 

0*2(0 



1 + au,-l{() + Ov,l c 2niC/e L + Q ^(C) 



o(l), and, by Corollary 5.1, one has = 0(t Vjll ). Therefore, to 



By (5.13), sup 
prove (5.31), it suffices to check that, for C G R 

(5.33) 5 (C) := = o(V). 

Let us check this. We know that 

(1) g is analytic in the half plane {ImC < Y} and tends to zero as ImC — ► — oo (as it is the 
sum of the Fourier series terms with the negative indexes of a function analytic in the strip 
{|ImC| <Y} ); 

(2) for Tm£ = Y, one has \g\ < C (by (5.13)). 

This implies that \g\ < Ce- 2 ^ Y ~ lraC -^ £ in the half plane {ImC < Y}. In view of (5.20) and (5.19), 
this implies (5.33), hence, (5.29). 
Finally, we check that 

(5.34) <p(r v ) = (p(b„,o) ~ <PM + 0(i 2 Kv t v , v , e~ 2nY / £ ). 

Therefore, for ( 6 1, we substitute the representations (5.21) and (5.31) in the relation b u = R u a u , 
and integrate the result over the period. As a v $ is the zeroth Fourier coefficient of a V) the mean value 
of the error term in (5.31) is zero. Hence, b Vt o = r u a U) o (l + 0{t 2 hv t %v , e- 2nY / e )) which implies (5.34). 
Representations (5.29), (5.34) and estimate (5.23) imply (5.22). The proof of Proposition 5.1 is 
complete. □ 

Proof of Lemma 5.1. We check the first and the second point for a u ; for b u , the proof is similar. 
Theorem 5.1 implies that, for |Im£| < Y, a u admits the representation 

(5.35) a u (() = a Ut0 (l+g )+a Utl e 2ni < :/£ (l+g 1 ) where |<? | + \gi\ = o(l). 

Therefore, the possible zeros of a u in the strip {|Im£| < Y} are located in o(e)-neighborhoods of the 
points 

(5.36) In (— a„ o/ov i) + /e, I € Z. 

This, Corollary 5.1 and the first point in Lemma 5.1 imply the second point of Lemma 5.1. 
To prove the first point of Lemma 5.1, we apply Rouche's Theorem to the functions / = ayfi+a u ^e 2m ^^ £ 
and 5f = a Ut ogo + a Vj \e 2m( *' £ g\. Clearly, all the zeros of / are simple and they are all listed in (5.36). 
Let C a be one of them. One compares / and 5f on the circle centered at £ a of radius ce (where c is a 
fixed positive sufficiently small constant independent of e). As 

/(C) l-n + l-l/ U ' 

then, on such a circle, one has 5f / f = o(l). This and Rouche's Theorem imply that a v has a unique 
simple zero in ce-neighborhood of Ca- This implies the first two points of Lemma 5.1 for a v . 

To prove the last point of Lemma 5.1, we compare the zeros of the functions b v b* v and a u a* inside the 
strip {— Y < ImC — 0}- We use the following observations: 
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• by the first two points of Lemma 5.1, in the strip {—Y < Im£ < 0}, all the zeros b u b* are 
zeros of b v , and all the zeros a v a* are zeros of a u ; 

• we know dyO^y — bj/b^ — det T u and that T u — 0(1) (see the second point of Theorem 5.1). 

So, the zeros of a v a* v have to be exponentially close to those of a v a* v — detT„, i.e. to the zeros of b v b* v . 
To study the distance between the zeros of a u a* and those of a u a* — detTj,, we again use Rouche's 
Theorem. Therefore, we pick £ a , a zero of a u and compare the functions / = a v a* u and 5f = detTj, on 
C r , the circle centered at Ca of radius 

r e 



where tq is a fixed positive constant, sufficiently large but independent of e. Note that, by Corollary 5.1, 
one has 

(5.37) \r\^r etl v (ReE). 

When applying Rouche's theorem, we have to control / on C r . Therefore, we use the relation 

(5.38) / / (0 = -^ ! <o^ ) o(l + o(l)) for |C - Ca| < e 2 - 

We prove (5.38) later, and, now, we use it to complete the proof of Lemma 5.1. By means of (5.38) 
and (5.37), for |£ — £ a | = r, we get 

2-7T 

1/(01 = — <o^,or(l + o(l)) = 27rr (l + o(l)). 



As 5f = detTjy = O(l), this implies that 

mo 



max 

IC-C«l=r 



< C/r . 



/(C) 

So, if ro is fixed sufficiently large, then, for sufficiently small e, / — 5f has one simple zero inside the 
circle |£ — ( a \ = r. As this is a zero of b v , and as r admits the estimate (5.37), this implies the third 
point of Lemma 5.1. 

To complete the proof of this lemma, we only have to check (5.38). Therefore, first, we note that, 
by (5.35), for — Y < Im£ < 0, one has 

a* u = a* (1 + o(l)) + a^e-^U + o(l)) = < (l + o(l) + o (|^)) = a* (1 + o(l)), 

a* 

where we have used Corollary 5.1 to estimate The result of this computation and (5.35) imply 
that, for -Y < lm( < 0, 

(5.39) /(C) = o* 0^,0(1 + o(l)) + a* a Vll ^ e (l + o(l)). 

The Cauchy estimates applied to o(l), the functions from (5.39) give J^o(l) = o(l) in any fixed 
compact of the strip {—Y < lm( < 0}. Therefore, for \( — ( a \ = e 2 , we get 

/'(C) = ^ a: i0 a,, ie 2 ^/ £ (l + o(l)) + o« a,, ) + o« a„ iie 2 ^). 

As a u ^e 2 ^/ £ = -nun, this implies (5.38). This completes the proof of Lemma 5.1. □ 

5.3. Asymptotics of the coefficients of the monodromy matrix. Using Theorem 5.1 and Propo- 
sition 5.1, we prove Theorem 2.2. Actually, we compute only the matrix M n corresponding to the 
consistent basis {f n ,f*}. The asymptotic representations for the coefficients of the matrix Mq are 
obtained similarly. The proof consists of two steps. 



5.3.1. Combinations of Fourier coefficients. First, we define the functions a u and the quantities 
T Vj u, Th, 6 and z u introduced in (2.15) - (2.18) in terms of the Fourier coefficients of the transition 
matrices. The asymptotics of the Fourier coefficient combinations met here are computed in terms of 
the iso-energy curve T in section 12. 

1. The phases. The phases <&„ are defined by the formulae 

= | (vCott.o) + f(ao,o) ~ p(bn,o) + ^( fo o,o)) , 
(5-40) . I 

®tt = 2 (^( a T,o) + ^(ao,o) + ~ <p(bo,o)) ■ 
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In section 12.2, we check that these phases admit the representations (2.19). They imply in particu- 
lar (2.28). 

2. The constant 9. Let 

(5.41) 9 = - ^z£ 

a-n-,0 

Note that, in view of (5.8), one has 



Qti-,0 

«o,o 



detTr 



o 



In section 12.1, we prove that 9 admits the representations (2.21) which, in particular imply (2.32). 
3. The coefficients Th and T v u . Let 



(5.42) 



Th = | Ott.o «o,o I 1 and T v 



a-u,o 



Using computations analogous to those done in section 12.1, one proves representations (2.20). These 
show that, for E G VI, 



(5.43) 

4- The constant z v . Let 
(5.44) 



\T h \ ^ \th\ = \th,oth,n\ and \\T V 



"2^ 



Using computations analogous to those performed in section 12.2, one proves (2.22). Estimate (2.23) 

is proved in the section 11.3. It implies (2.33). 

5. The functions a u . Define a u = a v ja v ^. One has 

Lemma 5.2. Fix < y < Y m . For sufficiently small e, in the case of Theorem 5.1, for |Im£| < y 
and E G V*, one has (2.18). 



Proof. Start with (5.32) or, equivalently, with 



(5.45) 



a 



v = l + g(0 + T v „e 2 ™^ (1 + g(0) , 9(0 



tv-i(C) 



, 9(0 



a »A0 



When proving (5.33), we have seen that \g\ < Ce~ 2 ^ Y - lmC -)h i n the half plane {Im( < Y}. Similarly, 
one proves that \g\ < Ce~ 2n( - Y+lm ^/ £ in the half plane {Im£ > — Y}. In view of (5.20) and (5.19), in 
the strip |Im£| < y, one has \T Vil/ e 2m< * £ \ < C. These three estimates imply (2.18). □ 
Note that (2.18) can be simplified into (2.26). 

6. Real analyticity. Note that T VtU , Th, 9 and z v , regarded as functions of E, are real analytic 
in V £ (this follows from the definitions of I • I and ¥?(•))• Therefore, each of them is invariant with 
respect to the operation * (see (2.10)). 

5.3.2. Computing the matrix M n . The representation (5.7) and the relation b v = R u a u imply that 
(5.46) A n = a n ao + R^R^a^aQ and = R§ a^ao + R n a n aQ. 

Now, for v G {0, 7r}, 



in (5.46), we substitute the representation a v = | a^o | e Vv "" / ' u -' a v ; 



• in (5.46), we replace the functions R v by their representations (5.17); 

• we express the Fourier coefficient combinations we meet in terms of <h u , T V)U , Th, 9 and z v ; 

• we use detToT^ = 1. 

• we use the invariance of <& u , T v>u , Th, 9 and z v with respect to the transformation *. 
This leads to the formulae 

o^e^Co , _ .».i*ir*a [9 + 1/9 , T h O w + O* , O w /9 + O* 9 O w O* 



(5.47) A„ = 2 ■ 
and 



+ a n Qfg e l ~ 



+ 



+ 



n 



(5.48) 



-tA 



2 — h a n e * 



a e l e 1 9 + a* Q e~ l e 9 a e* * O + a^e' 1 e O n 



+ 



n 
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where A = 99(60,0) — <£>(ao.o)- Furthermore, one checks that, for |Im£| < y and E G one has 

(5.49) a n a* {■■■}= U + ^ + O (pT V)0 , pT Vt7T , T h , pT Y /T h ) 
and 

(5.50) a n e iS ? {•■■} = \ ^ (I ^ + 9e"^ + O (pT vfi , pT Vj7T , T h , P T Y /T h ) . 

In (5.49) and (5.50), the terms with the curly brackets are the ones from (5.47) and (5.48) respectively, 
and p = p((/e). These two representations follow from estimates (2.26), (2.27), (2.28), (2.32), (5.18) 
and (5.43). We omit the elementary details. 

Finally, we "kill" the constant factor e~ tA in (5.48) by replacing the consistent basis {f^,/*} with 
the consistent base {g,g*} where g = e~ l ^/ 2 f n : for the monodromy matrix corresponding to {g,g*}, 
the coefficient with index 11 is equal to A w , and the coefficient with index 12 is equal to B n e~ lA . For 
the coefficients M\\ and M12 of this new monodromy matrix, we keep the old notations A n and B^. 
With this "correction", the asymptotic representation (2.15) follows from the representations (5.47) 
and (5.49), and the asymptotic representation (2.16) follows from the representations (5.48) and (5.50) 
This completes the proof of Theorem 2.2. □ 



6. Periodic Schrodinger operators 

In this section, we discuss the periodic Schrodinger operator (0.2) where V is a 1-periodic, real valued, 
L^ oc -function. First, we collect well known results needed in the present paper (see [13, 6, 19, 21, 26]). 
In the second part of the section, we introduce a meromorphic differential defined on the Riemann 
surface associated to the periodic operator. This object plays an important role for the adiabatic 
constructions (see [9]). 

6.1. Analytic theory of Bloch solutions. 

6.1.1. Bloch solutions. Let ip be a nontrivial solution of the equation 

d 2 

(6.1) — r? ip(x) + V(x)il)(x)=Si/}(x), xGR, 

ax z 

satisfying the relation tp {x + 1) = Xijj (x) for all x G R with A G C independent of x. Such a solution 
is called a Bloch solution, and the number A is called the Floquet multiplier. Let us discuss properties 
of Bloch solutions (see [13]). 

As in section 1.1, we denote the spectral bands of the periodic Schrodinger equation by [E±, E2], 
[E3, £4], . . . , [E2ra+i, E2n+2], ■ ■ ■ ■ Consider S±, two copies of the complex plane £ G C cut along the 
spectral bands. Paste them together to get a Riemann surface with square root branch points. We 
denote this Riemann surface by S. In the sequel, ir c : S 1— > C is the canonical projection. 

One can construct a Bloch solution ip(x,£) of equation (6.1) meromorphic on S. For any £, we 
normalize it by the condition ip(l,£) = 1. Then, the poles of £ 1— > ip{x,£) are projected by tt c either 
in the open spectral gaps or at their ends. More precisely, there is exactly one simple pole per open 
gap. The position of the pole is independent of x (see [13]). 

Let : : S 1— > S be the canonical transposition mapping: for any point £ G S, the point £ is the unique 
solution to the equation tt c (£) = E different from £ outside the branch points. 

The function x 1— ► ip(x, £) is one more Bloch solution of (6.1). Except at the edges of the spectrum (i.e. 
the branch points of S), the functions ip(-,£) and ip{-,£) are linearly independent solutions of (6.1). 
In the spectral gaps, they are real valued functions of x, and, on the spectral bands, they differ only 
by the complex conjugation (see [13]). 

6.1.2. The Bloch quasi-momentum. Consider the Bloch solution ip(x,£). The corresponding Floquet 
multiplier A (£) is analytic on S. Represent it in the form X(£) = exp(ik(£)). The function k(£) is 
the Bloch quasi-momentum. 

The Bloch quasi-momentum is an analytic multi-valued function of £. It has the same branch points 
as ip(x,£) (see [13]). 

Let D G C be a simply connected domain containing no branch point of the Bloch quasi-momentum k. 
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On D, fix ko, a continuous (hence, analytic) branch of k. All other branches of k that are continuous 
on D are then given by the formula 

fc±,i(£) = ±k (£) + 2ttZ, Z e Z. 

All the branch points of the Bloch quasi-momentum are of square root type: let E\ be a branch point, 
then, in a sufficiently small neighborhood of Ei, the quasi-momentum is analytic as a function of the 
variable \f£ — E\\ for any analytic branch of k, one has 

k{£) = h + ci^S -Ei + 0(8 - Ei), 
with constants k[ and q depending on the branch. 

Let C+ be the upper complex half-plane. There exists k p , an analytic branch of k that conformally 
maps C+ onto the quadrant {k £ C; Imk > 0, Re/c > 0} cut along compact vertical intervals, say 
ttI + Hi where I £ N* and C R, (see [13]). The branch k p is continuous up to the real line. It is 
real and increasing along the spectrum of Hq; it maps the spectral band [-E^n-i* -^2n] on the interval 
[ir(n— l),7rrt]. On the open gaps, Re k p is constant, and Imk p is positive and has exactly one maximum; 
this maximum is non degenerate. 

We call k p the main branch of the Bloch quasi-momentum. 

Finally, we note that the main branch can be analytically continued on the complex plane cut only 
along the spectral gaps of the periodic operator. 

6.2. Meromorphic differential Q. 

6.2.1. The definition and analytic properties. On the Riemann surface S, consider the function 

j}ijj(x,£) (ip(x,£ ) — ik(£ )xtp(x,£ ) ) dx 

(6.2) u(S) = \ ; '- . 

J ip(x,£)ip(x,£)dx 

where k is the Bloch quasi-momentum of ip, and the dot denotes the partial derivative with respect 
to E. This function was introduced in [10] (the definition given in that paper is equivalent to (6.2)). 
In [10], we have proved that to has the following properties: 

(1) the differential f2 = lj d£ is meromorphic on 5; its poles are the points of P U Q, where P is 
the set of poles of ijj(x,£), and Q is the set of points where k'(£) = 0; 

(2) all the poles of f2 are simple; 

(3) Vpe P\Q, respft = 1; Vg € Q\P, res g O = -1/2; VrGPnQ, res r ft = 1/2. 

(4) if 7T C (£ ) belongs to a gap, then lu(£) £ R; 

(5) if tt c (£) belongs to a band, then uj(£) = uj(£). 

The following quantities appeared in the description of the spectrum of H z e (see sections 1.6 and 1.7) 

(o.3) MV0= i( e „ (m _y, 

where 

(6.4) n (V) = exp(l n (V)), l n {V)= ! 0(f), 

J 9n 

and g n is a simple closed curve on S such that 

• g n is located on C\o~(Ho), the sheet of the Riemann surface S where the Bloch quasi-momentum 
of ip(x,£) is equal to k p (Tr c (£)) for lmir c (£ ) > 0; 

• ^c(gn) is a positively oriented loop going once around the n-th spectral gap of the periodic 
operator Hq. 

We prove 

Lemma 6.1. The integral l n is real valued. 

Proof. Let £q be a point that projects onto an internal point of a spectral band. Let U be a 
neighborhood of £q where is analytic. Let here £* = 7r~ 1 (7r c (i5)). By the fifth property of lu, for 
£ & U, one has uj(£) = to(£*). Consider g n , the integration contour for l n . We can and do assume that 
^cidn) (as a set, but not as an oriented curve) is symmetric with respect to the real line. As Tr c {g n ) 
intersects the real line at internal points of spectral bands, starting from one of these intersections, 
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one can continue £ \— ► £* to g n as a continuous function; then, one has uj(£) = w(£*). Note that 
5n = -ffn- One has 



(6.5) 



J n(£) = f u(£)d£ = ! u(£*)d£= f u(£)d£ 
•> an Jg n Jg'n Jg* n 

= - I uj{£) d£ = - I uj{£) d£ = - I Q(£). 

J On J On J On 



On S, there are exactly two points, say q and q, in Q that project inside the ra-th spectral gap of Hq. 
Furthermore, on S, there is exactly one point, say p, in P that projects inside the n-th spectral gap 
or at one of its edges. On S \ {q, q,p}, up to homotopy, one has 

g n = -9n+ ^2 

£ e{<j,g,p} 

where C{£) is a infinitesimally small, positively oriented circle centered at £. This and the description 
of the poles of Q imply that 

(6.6) I Q(£) = - f n(£) + 2m res £ n(£) = -f 

•Jgn Jgn ££{q,q,p} 9n 

Relations (6.5) and (6.6) imply that l n = l n . This completes the proof of Lemma 6.1. □ 
Lemma 6.1 imply 

Corollary 6.1. One has 6 n (V) > and A n (V) > 1. 



7. The consistent solutions 

In this section, we describe the consistent solutions (fv)ve{o,Tr} used in section 5.1. Many of the results 
presented here are taken from [9]. 

In [11] and [9], we have developed a new asymptotic method to study solutions to an adiabatically 
perturbed periodic Schrodinger equation i.e., to study solutions of the equation 

(7.1) C) + (F(x) + W{ - £X + C))V,(x ' C) = E ^ C) 

in the limit e — > 0. The function £ i— > W(Q is an analytic function that is not necessarily periodic. 
The main idea of the method is to get the information on the behavior of the solutions in x from 
the study of their behavior on the complex plane of £. The natural condition allowing to relate the 
behavior in £ to the behavior in x is the consistency condition (5.3): one can construct solutions 
to (7.1) satisfying this condition and having simple standard asymptotic behavior on certain domains 
of the complex plane of £. 

We first describe the standard asymptotic behavior of the solutions studied in the framework of the 
complex WKB method. The domains where these solutions have the standard behavior are described 
in terms of Stokes lines. So, next, we describe the Stokes lines for V, W and E considered in this paper. 
Finally, we describe fo and /„■, the solutions used to construct the consistent bases and transitions 
matrices of Theorem 5.1. 



7.1. Standard behavior of consistent solutions. We now discuss two analytic objects central to 
the complex WKB method, the complex momentum defined in (1.1) and the canonical Bloch solutions 
defined below. For £ G T>(W), the domain of analyticity of the function W, we define 

(7.2) £(() = E-W(() 

The complex momentum and the canonical Bloch solutions are the Bloch quasi-momentum and par- 
ticular Bloch solutions of the equation 

(7.3) _^ + ^ = f( C )^. 



considered as functions of £. 
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7.1.1. The complex momentum. For £ G 2?(W), the domain of analyticity of the function W, the 
complex momentum is given by the formula = k (£(C)) where A; is the Bloch quasi-momentum 
of (0.2). Clearly, k is a multi- valued analytic function; a point C such that W'(C) ^ is a branch 
point of k if and only if 

(7.4) Ej + W(C) = E for some j G N*. 
All the branch points of the complex momentum are of square root type. 

A simply connected set D C T>(W) containing no branch points of k is called regular. Let k p be a 
branch of the complex momentum analytic in a regular domain D. All the other branches analytic in 
D are described by 

(7.5) K^j = ±k p + 27rm where m G Z. 

7.1.2. Canonical Bloch solutions. To describe the asymptotic formulae of the complex WKB method, 
one needs Bloch solutions of equation (7.3) analytic in ( on a given regular domain. We build them 
using the 1-form Vt = uj d£ introduced in section 6.2. 

Pick (q, a regular point. Let £q = 5 (Co)- Assume that £q G" P U Q (the sets P and Q are defined 
in section 6.2). In Uq, a sufficiently small neighborhood of £q, we fix k, a branch of the Bloch 
quasi-momentum, and ip±(x,£), two branches of the Bloch solution i/j(x,£) such that k is the Bloch 
quasi-momentum of Also, in Uq, consider £l±, the two corresponding branches of f2, and fix 

a branch of the function £ i— ► q{£) = ^k'(£). Assume that Vq is a neighborhood of Co such that 
£(V ) C Uq. For C G Vb, we let 

(7.6) y±(x,() =q(£)e I ^ n± 'ip ± (x,£), where £ = £(C)- 

The functions ^± are the canonical Bloch solutions normalized at the point Co- Its quasi-momentum 
is K (() = k(E-W(0). 

The properties of the differential O imply that the solutions *5/± can be analytically continued from 
Vq to any regular domain D containing Vq. 
One has (see [11]) 

(7.7) w^+(; c), *_(-, 0) = «;(*+(-, Co), *-(-, Co)) = fc / (£oMV'+(-,£o),^-(-,£o)) 

As £ Q U {E 1 ;, Z > 1}, the Wronskian «;(*+(■, C), 0) does not vanish. 

7.1.3. Solutions having standard asymptotic behavior. Fix E 1 = £?o- Let D be a regular domain. Fix 
Co G so that £(Co) -P U Q. Let k be a branch of the complex momentum continuous in D, and 
let ty± be the canonical Bloch solutions defined on D, normalized at Co an d indexed so that k be the 
quasi-momentum for 

We recall the following basic definition from [9] 

Definition 7.1. Fix s G {+,—}• We say that /, a solution of (7.1), has standard behavior (or 
standard asymptotics) / ~ exp(s^ ndQ) ■ ^ s in D if 

• there exists Vq, a complex neighborhood of Eq, and X > such that / is defined and satis- 
fies (7.1) and (5.3) for any (x,(,E) G [-X,X] x D x V ; 

• / is analytic in C G D and in E G Vq; 

• for any K, a compact subset of D, there is 1/ C Vo, a neighborhood of Eq, such that, for 
(x,£,E) G [— X, X] x K x V , f has the uniform asymptotic 

f = e s U CKd ((y s + (l)), as e^O, 

• this asymptotic can be differentiated once in x without loosing its uniformity properties. 

db- f' ndC 

Let (/+,/_) be two solutions of (7.1) having standard behavior /± ~ e e J ' in D. One 
computes 

«;(/+, /_) = tu(*+,*_) + o(l). 

By (7.7), for C in any fixed compact subset of D and e sufficiently small, the solutions (/+,/_) are 
linearly independent. 
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7.2. Complex momentum and Stokes lines for W(C) = acos£. Now, following [9], we discuss 
the complex momentum and describe the Stokes lines for V, W and E considered in this paper. In 
particular, from now on, we assume that 



(7. 



W(C) = acos(C) hence, £(£) = E — acos(C), 



that E belongs to J, a compact interval satisfying the condition (TIBM) from section 1.2, and that 
all the gaps of the periodic operator Hq are open. 

7.2.1. Complex momentum. 1. The branch points of the complex momentum are located on the 
lines of the set arccos (R) which consists of the real line and the lines {Re £ = ttI} for I £ Z. The 
set of branch points of k is 2-/r-periodic and symmetric with respect both to the real line and to the 
imaginary axis. 

Dehne the half-strip II = G C; < Re£ < ir, Im£ > 0}. It is a regular domain. Consider the 
branch points located on dll, the boundary of II. £ bijectively maps 911 onto the real line. So, for 
any j £ N, there is exactly one branch point solution to (7.4); we denote it by Under condition 
(TIBM), the branch points Qm and C2n+i are located on the interval (0,7r), i.e. < ^2n < C2n+i < tt- 
The branch points (j, C2, ••• (,2n-i are located on the imaginary axis and satisfy < Im £271-1 < 
••• < Im^2 < I m Ci- The other branch points are located on the line {Re£ = 7r}, and one has 

< Im C2n+2 < I™ C2n+3 < In Fig. 7, we show some of these branch points. 

2. £ conformally maps the half-strip II onto the upper half 
of the complex plane. So, on II, we can define a branch of the 
complex momentum by the formula 



(7.9) 



K p (ip) = k p (E — a cost/?), 



k p being the main branch of the Bloch quasi-momentum for 
the periodic operator (0.2). We call k p the main branch of the 
complex momentum. 

The discussion in section 6.1.2 implies the following. First, k p 
conformally maps II into the first quadrant of the complex 
plane. Fix I, a positive integer. The closed segment Z[ := 
[C2Z-liC2i] C <9II is bijectively mapped on the interval [ir(l — 
1) , 7rZ] ; on the open segment gi := ((21X21+1) C dll, the real 
part of k equals to irl, and its imaginary part is positive. Two 
of the intervals (zi)i and (gi)i are shown in Fig. 7. 

7.2.2. Stokes lines. Let (q be a branch point of the complex 
momentum. 

A Stokes line beginning at £0 is a curve 7 defined by the 
equation Im J^ (k (£) — K (Co))d£ = (where k is a branch 
of the complex momentum continuous on 7). There are 
three Stokes lines beginning at each branch point of the 
complex momentum. The angles between them at the 
branch point are all equal to 2n/3. 

Let us discuss the set of Stokes lines for W(() = acos(. 
Due to the symmetry properties of £, the set of the Stokes 
lines is 27r-periodic and symmetric with respect to both 
the real and the imaginary axes. So, it suffices to describe 
the Stokes lines in II. Here, we follow [9]. 

In Fig. 8, we have represented Stokes lines in II by dashed 
lines. 
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Figure 8: The Stokes lines 



Elementary properties of Stokes lines. Recall that the 

ends of the intervals (gi)i are branch points and, reciprocally, any branch point located on 911 is an 
end of one of the gis. 

Consider the Stokes lines beginning at the ends of g n . The right end of g n is C2n+i • One of the Stokes 
lines beginning at this point goes to the right of C2n+i along K; the two other Stokes lines beginning 
at C2n+i are symmetric with respect to the real line. Similarly, one of the Stokes lines beginning at 
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(,2m the left end of g n , goes to the left of ^2n along R, and the two other Stokes lines beginning at £271 
are symmetric with respect to the real line. 

Consider the Stokes lines beginning at the ends of g\ for either l>n-\-loxl<n — 1. One of these 
Stokes lines coincides with g\. Let Co be one of the ends of g\. The two Stokes lines beginning at Co 
and different from g\ are symmetric with respect to the line {Re£ = Re Co}, see Fig. 8. 

Global properties of the Stokes lines in IT. First, we discuss the Stokes lines starting at C2n+i,---, C2n+4 
and (,2n denoted respectively by "a",..., "d" and "e". They are shown in Fig. 8 and described by 

Lemma 7.1 ([9]). The Stokes lines "a",..., "d" and "e" have the following properties: 

• the Stokes lines "a" and "e" stay inside LT, are vertical and disjoint; 

• the Stokes line "c" stays between "a" and the line {Re£ = tt} (without intersecting them) and 
is vertical; 

• before leaving LT, the Stokes lines "b" stays vertical and intersects "a" at a point with positive 
imaginary part; 

• before leaving II, the Stokes lines "d" stays vertical and intersects "c" above C2n+3, the beginning 
of "c". 

The term "vertical line" used in this lemma means a smooth curve intersecting the lines {Im£ = C} 
transversally. The proof of Lemma 7.1 can be found in [9]. 

Now, consider the Stokes lines located in LT and starting at C2n-i, C2n-2 and C2n.-3- We respectively 
denote them by "f" , "g" and "h" , see Fig. 8. One proves 

Lemma 7.2. The Stokes lines "f", "g" and "h" have the following properties: 

• the Stokes line "g" is vertical and stays between "e" and the line {Re£ = 0} without intersecting 
them; 

• before leaving LT, the Stokes lines "f" stays vertical and intersects "e" at a point with positive 
imaginary part; 

• before leaving IT, the Stokes lines "h" stays vertical and intersects "g" above Qm-2, the beginning 
of "g". 

We omit the proof of this lemma as it is similar to the proof of Lemma 7.1. 

7.3. Two consistent solutions. We now introduce two solutions of (5.2) satisfying (5.3). For y > 0, 
we define S y = {|Im£| < y}. 

Fix Y > Im£2n+4- The solutions we describe are analytic in the strip Sy- 

We first describe the branch of the complex momentum used to write the asymptotics of these solutions. 
Define the strip 

S p = {( G C;0 < ImC < min(ImC2„-i,ImC2n+ 2 )}. 

It is regular. Analytically continue k p to this strip. Recall that the integer n in the condition (TIBM) 
is even. Let 

(7.10) k(C) = MO - nir. 

As n is even, the discussion in the section 7.1.1 shows that k is a branch of the complex momentum. 
It is continuous up to the boundary of the strip S p ; one has 

K((2n) = «(C2n+l) = 0. 

7.3.1. The solution f w . Consider V n , the subdomain of the domain D n = {|Im (\ < Y, < Re ( < 2n} 
shown in Fig. 9(a). Its boundary consists of the lines bounding and of the segments of Stokes lines 
and lines {Re£ = C} beginning at the intersection points of Stokes lines. The domain T>^ is simply 
connected. 

Let K n be the analytic continuations of k from S p to T> n , i.e., for £ € T> n n S p 

(7.11) k*(C) = k(0- 

(tt) 

Let \&Y be the canonical Bloch solution analytic T> n , normalized at tt and such that n n is its Bloch 
quasi-momentum. In [9], we have proved 
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Proposition 7.1 ([9]). Fix E = E* £ J. For sufficiently small e, there exists f n , a solution to (5.2), 
satisfying (5.3) and analytic in the strip S Y that, on T> n , has the standard behavior 



fC 

(7.12) U ~ exp ( - / K w d( ) 



7.3.2. TTie solution Jq. Consider 2?0) the subdomain of the domain Dq = {|Im£| < Y, —ir < Re£ < 7r} 
shown in Fig. 9(b). Its boundary consists of the lines bounding Dq and of the segments of Stokes lines 
and lines {Re£ = C} beginning at the intersection points of Stokes lines. The domain T>§ is simply 
connected. 

Let kq be the analytic continuation of — k from S p to T>q i.e., for £ € T>o H S p , 
(7.13) «o(C) = -k(C). 

Let be the canonical Bloch solution analytic Vq, normalized at and such that kq is its Bloch 
quasi-momentum. One has 

Proposition 7.2. Fix E = E* G J . For sufficiently small e, there exists fo, a solution to (5.2), 
satisfying (5.3) and analytic in the strip S Y that, on T>q, has the standard behavior 

i f C -A -(o) 



(7.14) f ~ exp I - K d( 1 * 

The proof of Proposition 7.2 is similar to that of Proposition 7.1; we omit it. 



8. TWO CONSISTENT BASES 

In this section, we construct the consistent bases used in the section 5.1. 

Fix v £ {0,7r}. The solution /* is related to f u by the transformation (2.10). First, we compute the 
asymptotics of /*. Then, we compute the asymptotic of the Wronskian w(f u ,f*). This Wronskian is 
constant up to a factor of the form (1 + o(l)). Finally, we correct / so that 

(1) the Wronskian w(f u ,f*) be constant (and, thus, /„ and /* form a consistent basis), 

(2) the "new" solutions /„ and /* have the "old" behavior in the strip Sy- 

The constructions described here are standard for the adiabatic complex WKB method. The proofs 
of Lemmas 8.1 and 8.2, and of Theorem 8.1 below essentially repeat the proofs of the analogous 
statements found in [10] and are therefore omitted. 

8.1. Asymptotics of /*. To discuss the asymptotic behavior of /*, we need some additional material. 

8.1.1. Preparation. Define 30 = ( — C2m C2n) C M. and 3^ = (C2n+i> 27r — C2n+i) C R. Note that 8 maps 
30 into the ra-th spectral band, and 3^ into the (n + l)-st spectral band. 

Recall that the leading terms of the asymptotics of the solutions having standard asymptotic behavior 
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are fixed up to the choice of q, the branch of ^/ A;' (£(£)) from the definition of the canonical Bloch 
solution, see (7.6). Let q v be the branch of q from the definition of . Fix it so that 

q u (C) > for C e l v . 

This choice is possible as, inside any spectral band of the periodic operator, k' p > 0. 

8.1.2. The asymptotics. Let T>* be the domain symmetric to T> v with respect to the real line. Note 
that j.cD.nD*. One has 

Lemma 8.1. In T>*, the solution f* has the standard behavior 

(8.1) /*~e~i/^* d ^ ) '*(x,C)- 

Here, is the branch of the complex momentum analytic in T>* that coincides with k v on $ u ; the 
function v|/v'* is the canonical Bloch solution analytic in T>* that coincides with (complementary 
to from the asymptotics of f u ) on % v . 

The proof of Lemma 8.1 mimics that of Lemma 6.1 in [10]. 
Note that = <, and that *^ } '* = (*f } )*- 

8.2. The Wronskian of f v and /*. The solution f v and /* are analytic in the strip S y ; so does 
their Wronskian. As both f v and /* satisfy condition (5.3), it is e-periodic in £. One has 

Lemma 8.2. The Wronskian of f u and f* admits the asymptotic representation: 

(8.2) w(f u JZ) = w(*M,vM)\ <=t , + gu, ( € S y . 

Here, g v is a function analytic in S y such that, for real £ and E, Reg^, = 0. Moreover, g v = o(l) 
locally uniformly in any compact of S y provided that E is in a sufficiently small complex neighborhood 
ofE . 

The proof of Lemma 8.2 mimics that of Lemma 6.2 in [10]. 

Remark 8.1. Note that \ ^ ) )\ c=u / as 8(u) PUQ. 

As g v , the error term in (8.2), may depend on £, we redefine the solution f v setting 



f»:=fv/Q where Q = yll + g/w^M-) 



In terms of this new solution /„, we define the new /*. The solutions (fuifu) form the basis the 
monodromy matrix of which we study. For these "new" f v and /*, we have 

Theorem 8.1. The solutions f v and f* are analytic in S y , satisfy the condition (5.3), and 

(8.3) w{uj:)=w{*¥ ) M-)k=»- 

Moreover, f v has the standard behavior, (7.14) or (7.12), inT> v , and f* has the standard behavior (8.1) 
in V*. 

The proof of Theorem 8.1 mimics that of Theorem 6.1 from [10]. 

Let £ i — > ip±(x, £(0) be the two branches of the Bloch solution £ i— ► ip(x,£(Q)) that are analytic 
in Q £ S p and such that k, the branch of the complex momentum defined in the beginning of the 
section 7.3, is the Bloch quasi-momentum for By (7.7) and the definitions of the canonical Bloch 
solutions fy±, one has 



.4) w{^ , ^)\ c=u = s{v)k' p {£{u))w(i> + ^-)\ Q=v , where s{v) 



■1 if u = 0. 



9. Transition matrices 



In this section, we compute the asymptotics of the transition matrices T u defined by (5.6) for the bases 
(/„,#) for i/G{0,7r}. 
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9.1. Elementary properties of the transition matrices. One can easily express the coefficients of 
the transition matrices, see (5.6), via the Wronskians of the basis solutions; formulas (5.6) immediately 
imply 

Lemma 9.1. One has 

, on „ r ^_ ^(-.C + 27r),/ *(-,O) , , n _ w(fo(;QJA;C + 2*)) 

[ j 7rlU M/o(-,0,/o*(-,0) ' 7riU w(/o(-,C),fi(-,0) ■ 

and 

, Q - rn ^(/o(-,o,/:(-,o) , rn ^(/ 7 r(-,o,/o(-,o) 

M ao(C) ~w(U, OJ*MY boiC) - W (U,C),ttU)Y 

For v G {0, 7r}, by the definition of the standard behavior, the basis {/^, /*} is defined and analytic 
for (£, E 1 ) G x V(Y) where V(Y) is a neighborhood of E* G J; this neighborhood is independent 
of e. One has 

Lemma 9.2. Pick v G {0, tt}. The matrix T v is analytic and e-periodic in Q G Sy and analytic in 
E G V(Y). Moreover, detT u is independent of C, and does not vanish. 

Proof. As the solutions f v and /* are analytic functions of the variables C, and E, so are the Wronskians 
in formulae (9.1) and (9.2). Moreover, by (8.3), the Wronskians in the denominators of (9.1) and (9.2) 
do not vanish. This implies the analyticity of the coefficients of the transition matrices. The periodicity 
in £ follows from the fact that all the solutions satisfy (5.3). Finally, relations (5.6) imply that 

(9.3) w(U(x, C + 2vr), f^x, ( + 2vr)) = det T n w(f (x, C), ft{x, ())■ 

Now, (8.3) imply that detT n is independent of £. Similarly one checks that det To is independent of 
This completes the proof of Lemma 9.2. □ 

9.2. The asymptotics of the transition matrices. We first introduce some notations: 

(1) For the Fourier coefficients of a v and b v we use the notations introduced in (5.9), and recall 
that p(z) = e 2n \ Jmz \. 

(2) Let Y n , Y v ^ and Yq, Y V) q be the distances marked in Fig. 9(a) and Fig. 9(b) respectively. 
E.g., Yq is the imaginary part of the point of intersection of the Stokes lines "g" and "h" (see 
Lemma 7.2). Note that, for any u G {0,7r}, one has 

(9.4) < Y VtV <Y U <Y. 

(3) We use the branch k introduced in the beginning of the section 7.3; ip± (resp. £l±) are the 
branches of ip(x,S(-)) (resp. 0) such that k is the Bloch quasi-momentum of ip+. When 
integrating n (resp. integrating Q or continuing analytically ip) along a curve, we choose a 
branch of k (resp. f2, ip) near the starting point of the curve and then continue it along the 
curve. 

(4) Let 7 be a curve and g be a function continuous on 7. We denote by Aargo| 7 the increment 
of the argument q(Q along the curve 7. 

The asymptotics of the transition matrices coefficients are described by 

Proposition 9.1. Pick v G {0, tt}. Fix Y so that Y VtU < Y < Y v . There exists V V {Y), a complex 
neighborhood of E* independent of e, such that, for e sufficiently small, j G {0,1} and E G V U (Y), 
one has the uniform asymptotics 

( i f' 2tt(tt — v)i f \ 

(9.5) a l/ j=expls- / ndQ — j h / il s + «A arg q\ a + o(l) j , a = a u j, 



/if 2tt(tt — u)i f \ 

(9.6) b vj = exp ys- J^KdC- j^^~ e ^ + y^fi s + iAargg| / 3 + o(l)J , (5 = f3 v 

where s = +1 if v = tt, and s = — 1 if v = 0. 

In (9.5) and (9.6), one integrates along the curves shown in Fig. 10 chosen such that £(Q (-PUQ) 
along them; for each of the integration curves, q denotes a branch of ( 1— > \/k'(£ (£)) continuous on 
this curve. 

For (C,E) G S(Y) x V U (Y), one has the uniform estimates 

(9.7) av,-i{.0 = °( a v,o), = o(b u , ), a Vj2 (0 = o(p((/e)a Vil ), b v>2 (0 = o{p{C,/e)b u ^). 
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Figure 10: The integrations paths for Theorem 5.1 



In the remaining part of the present section, we first explain how Theorem 5.1 is deduced from 
Proposition 9.1. Then, we turn to the proof of Proposition 9.1. We begin with describing general 
asymptotic formulae for the Wronskians of two solutions having standard behavior; this material 
mostly stems from [10]. Then, using these formulae, we compute the Wronskians in the formulae for 
the transition matrix coefficients (see Lemma 9.1) and, thus, complete the proof of Proposition 9.1. 
Note that we carry out the analysis only for the asymptotics and the estimates for ao and bo. The 
coefficients a n and b n are analyzed in a similar way. 

9.3. The proof of Theorem 5.1. In section 10, we study the actions {S V) u)ve{o^} and prove 

Lemma 9.3. Pick E* G J. For each v S {0, tt}, one has S V)U (E*) = 2ttY v ^ u {E*). 

Lemma 9.3 and the condition (T), see section 1.5, imply that, in Proposition 9.1, we can choose Y so 
that (1) 2ttY > maxBgj Sh(E) and (2) Y v>v < Y < Y v simultaneously for v = and v = tt. We then 
define V* = Vq n V n . With this, each of the basis solutions /o, /q , f n and /* is defined and analytic in 
the domain (5.5). This and Lemma 9.2 imply the first and the second point of Theorem 5.1. 
In section 11, we derive the estimates of the third point of Theorem 5.1 from the asymptotics (9.5) 
and (9.6). 

Finally, the last point of Theorem 5.1 is an immediate consequence of the estimates (9.7). So, Theo- 
rem 5.1 is proved. 

9.4. General asymptotic formulae. We recall results from section 8 of [10]. Consider equation (7.1) 
assuming only that W is analytic and that E is fixed, say E = E$. Let h and g be two solutions 
of (7.1) having the standard asymptotic behavior in regular domains Dh and D g : 

(9.8) h~e< S i Kh ** h (x,Q, g~e'^ K3d(; y g (x,0. 

Here, Kh (resp. K g ) is a branch of the complex momentum analytic in Dh (resp. D g ), (resp. *S> g ) 
is the canonical Bloch solution defined on Dh (resp. D g ) and having the quasi-momentum Kh (resp. 
K, g ), and Qh (resp. Q g ) is the normalization point for h (resp. g). 

As the solutions h and g satisfy the consistency condition, their Wronskian is e-periodic in £. First, 
following [10], we describe the asymptotics of this Wronskian and of its Fourier coefficients. Then, we 
develop simple tools to compute some constants coming up in these formulae. 

9.4.1. Asymptotics of the Wronskian. Let d be a simply connected domain such that d C Dh n D g . 
Arcs. Let 7 be a curve connecting Q g to Qh going first from Q g to some point in d while in D g and, 
then, from this point to Qh while in Dh- We call 7 an arc associated to the triple (g,h,d) and denote 
it by 7(5, h,d). 

Two arcs associated to one and the same triple are called equivalent. 
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Continue Kh and n g analytically along 7(5, h,d). Then, there exists m G Z and a G {— 1,+1} such 
that for £ close to 7, one has 

(9.9) K fl (C) =<™ h (C)+27rm. 

We call cr = o~(g, h, d) the signature of 7, and m = m(g, h, d) the index of 7. These two integers do not 
change when we replace the arc 7 by an equivalent one. 

Meeting domains. A domain d is called a meeting domain if the functions IniK/j and Im.K g do not 
vanish and are of opposite sign in d. One has 

Lemma 9.4 ([10]). Suppose the functions Im/t/, and Iva.K g do not vanish in d. Then, d is a meeting 
domain if and only if a(g, h,d) = — 1. 

Fourier coefficients. Let S(d) be the smallest strip of the form {C\ < Im£ < C2} containing the 
domain d. One has 

Proposition 9.2 ([10]). Fix Eq. Let d = d(h, g) be a meeting domain for h and g, and m = m(g, h, d) 
be the corresponding index (at energy Eq). Then, there exists Vq a neighborhood of Eq such that for e 
sufficiently small, for E G Vq and C G S(d), the Wronskian of h and g is given by the formulae 

(9.10) w(h, g) = w m e^K-O^i + (i)), 
where 

(9.11) w m = (q g /q h )\ c=c , expf- f K g d( + [ Sl g ] w{^> + ,^. 

\ £ J 7(9 ,M) J( g J 

In these formulae: 

• w m is independent of C; 

• we choose the arc ^y(g,h,d) so that, along it, £(C) PUQ; 

• C l— > Qg(C) = y/k'(£(0) an d C l— * ^a(C) = Qg(£(0) are ^ e analytic continuations of the 
function and the 1-form from the definition of ' \& g along j(g,h,d). 

• ^ + = and is the canonical Bloch solution "complementary" to 

Fix K, a compact subset of S(d) . Then, there exists V a neighborhood of Eq in Vq such that the 
asymptotics (9.10) is uniform in K x V . 

The factor w m is the leading term of the asymptotics of the m-th Fourier coefficient of w(h,g). 

9.4.2. Closed curves and the index m. In practice, it is not too difficult to compute the index m. 
However, as one needs to control several Fourier coefficients of each Wronskian, the computations 
become lengthy. Fortunately, there is an effective way to compare the indexes of two (non-equivalent) 
arcs. To this end, we define the index of a closed curve. 

Closed curves. Let c be an oriented closed curve containing no branch points of the complex mo- 
mentum. Pick Co G c. In Vo, a regular neighborhood of Co, fix k, an analytic branch of the complex 
momentum. We call the triple (c, £o> k) a loop. 

We shall consider c as disjoint at Co and speak about its beginning and its end. Continue k analytically 
along c. This yields a new branch of the complex momentum in Vo- Denote it by k\ c . Hence, there 
exists a G { — 1, +1} and m G Z such that, for C G Vo 

(9.12) k| c (C) =<™(0 +2vrm. 

The numbers a = a(c, Coj k ) and m = m(c, Co; K ) are called the signature and the index of the loop 
(c,Co,«)- 

Consider two loops (ci,Ci,^i) and ^2X2,^2)- Assume that one can continuously deform c\ into C2 
without intersecting any branching point. Assume moreover that, in result of the same deformation, 
Ci becomes £2- This deformation defines an analytic continuation of K\ to a neighborhood of C%- If 
this analytic continuation coincides with K2 , we say that the loops are equivalent. The indexes m and 
a calculated for equivalent loops coincide. 

Let us explain how to compute the indexes m and a. Let G be the pre-image with respect to £ of the 
set of the spectral gaps of the periodic operator (0.2). Note that 

• on any connected component of G, the value of the real part of the complex momentum is 
constant and belongs to {71^; I G Z}; 

• locally, outside {C; W'(C) = 0}, all the connected components of G are analytic curves. 
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Now, we can formulate 



Lemma 9.5. Assume that c does not start at a point of G. Assume moreover that c intersects G 
exactly N times (N < oo) and that, at the intersection points, W ^ 0. Let r±, r2, ■ ■ ■ , r/v be the 
values that Re k takes consecutively at these intersection points as C moves along c (from the beginning 
to the end). Then, 

(9.13) cr(c, Co, «) = (-1) N , and m(c, Co,«) = - (tn ~ r N -i +r N -2 h (-l^Vi) . 

7T 

The proof of Lemma 9.5 mimics the proof of Lemma 8.2 in [10] which describes the index of a 2tt- 
periodic curve. 

Comparing the indexes of arcs. Let d and d be two (distinct) meeting domains for the solutions h and 
g, and let 7 and 7 be the corresponding arcs. One can write 

(9.14) 7 = c + 7) 

where c is a closed regular curve; its orientation is induced by those of 7 and 7. 

As a(g,h,d) = a(g,h,d) = —1, one has a{c,C g ,K g ) = 1. As an immediate consequence of the 
definitions, we also get 

(9.15) 771(3, h, d) = m(c, C g , K g ) + m(g, h, d). 

This formula and Lemma 9.5 give an effective way to compute the indexes of arcs. 

9.5. The asymptotics of the coefficient bo- The coefficient 60 of the matrix To is given in (9.2). 
As w(f n ,f*) is given by formula (8.3), we have only to compute w(f w (-, (), /o(-, C))- O ne applies the 
constructions of section 9.4 with 

h(x,() = /ttOcjO) g(x,0 = fo(x,(), D h = T> n , D g = V ; 

(9.16) Ch = C 9 = 0; 

(9.17) K h {() = k(C) for C ~ vr, and « ff (C) = -k(C) for ( ~ 0. 

In (9.17), k is the branch of the complex momentum defined in (7.10). 
Let Yq and Y v $ be the distances marked in Fig. 9(b). They satisfy (9.4). 

9.5.1. The asymptotics in the strip {—1^,0 < Im£ < Yo}- Let us describe do, the meeting domain, 
and 7(/o, fn, do), the arc used to compute w(f n , fo) in the strip 

s = {cec; -y„, <imc<r }- 

The meeting domain do- It is the sub domain of the strip So between the lines 71 and 72 defined by 

• the line 71 consists of the following lines: the Stokes line "e" symmetric to the Stokes line "e" 
with respect to the real line, the segment [0, C2n] of the real line, the segment [0, C2n-2] of the 
imaginary axis and the Stokes line "g" (see Fig. 8); 

• the line 72 consists of the following lines: the Stokes line "a" symmetric to the Stokes line "a" 
with respect to the real line, the segment [^n+ij^r] of the real line, the segment [7r,C2n+3] of 
the line Re£ = ir and the Stokes line "c" (see Fig. 8). 

The Stokes lines mentioned here are described by Lemmas 7.1 and 7.2. In particular, these lemmas 
imply that 71 n 72 = 0. 

Note that do does not intersect Z, the pre-image of the set of the bands of the periodic operator (0.2) 
with respect to the mapping £. So, in do, one has Ims 7^ 0. 

The arc 7(3, h,do). It is the curve /3o,o shown in Fig. 10; it stays in do and connects Cg = to Qh = tt- 
Index m. In view of (9.17), in do, one has Kh = —K g . This implies that m(g,h,do) = 0. 
The result. Proposition 9.2, formulae (9.2) and (8.3) imply that, for £ G So, 

(9.18) 6o = 6o(l + o(l)), 6 = exp^-^ J Kd( + J 0_ + iA arg q\^j where (3 = /3 ,o, 

and, as q, one can take any branch of the function ( 1— ► ^/k 1 {£((,)) continuous on (3. 
When deriving the formula for bo, we have used the facts that 

• £l g is the branch of f2_ corresponding to the branch k chosen above; 

• (Qg/Qh)(Ch) = e iAarg ^ 1 ' 3 as, at Ch, qh is real and \q g /q h \ = 1; 
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• the quantity Aargg^ does not depend on the branch of q 9 as long as it is continuous on /3. 

9.5.2. The asymptotics in the strip {—Yq < Im£ < — Y v ^}. Let us describe d\, the meeting domain, 
and 7 (/,[■, /o, di), the arc used to compute w(f n , fo) in this strip 

5i={(GC; -Y <Im(<-Y vfi }. 

The meeting domain d\ . Let d\ be the subdomain of the strip Si located between the Stokes line "a" 
(symmetric to "a" with respect to the real line) and 73, the curve which consists of the following lines: 

• the Stokes line "f" symmetric to the Stokes line "f" with respect to the real line, the segment 
[C2n-ii C,2n-2\ of the imaginary axis, and the Stokes line "g" symmetric to the Stokes line "g" 
with respect to the real line. 

The domain d\ is a meeting domain in view of 
Lemma 9.6. In d\, one has Im.K n = —Iuikq > 0. 

Proof. The sign of Ims remains the same in any regular domain D such that DflZ = fl. Moreover, 
the sign of Im k flips as £ intersects (transversally) a connected component of Z at a point where W 
does not vanish. 

By (7.10) and (7.11), one has Iuik^ = Ihik; = Iid.k p > in D n n IL As one goes from II to d\ in 
Ptt without intersecting Z, we get Im /«„-(£ ) > for Q G d\. Similarly, by (7.10) and (7.13), one has 
IniKo = -IniKp < in T>q D IT. Furthermore, to go from II to d\ staying in T>q, one has to intersect 
two connected components of Z, namely, the segment [— (,2niC2n] °f the real line and the segment 
[ — C2n-li C2n-i] of the imaginary axis. Hence, ImKo(C) < f° r C G d\. This completes the proof of 
Lemma 9.6. □ 

The arc 7(5, h,d\). It is the curve /?o,i shown in Fig. 10; it connects Cg = to Ch = tt- 

Index m. One has 

l(g,h,di) = co + 7(5, h, do), 
where cq is the closed curve shown in Fig. 11. By (9.15), we get 

m(g,h,di) = m(c ,0, k 9 ) + m(g,h,d ) = m(c ,0,n g ). 



2n+2 



So, the index m(g, h, d±) is equal to the index of the loop 
(co, 0,K g ). Recall that the indexes of equivalent loops co- 
incide. To compute the index, we pick a point (0 £ cq 
as shown in Fig. 11 and we replace the loop (cq, 0, k 9 ) by 
the equivalent loop defined by the same curve Co and the 
point Co- The branch of the complex momentum fixed 
for this new loop is the analytic continuation of the old 
branch along Co from to £0 in the clockwise direction. 
For this new branch, we keep the old notation n g . 
In view of Lemma 9.5, it is sufficient to compute k 9 at 
the intersections of cq and G. The set G is 27r-periodic 
and symmetric with respect to the real line and to the 
imaginary axis. The connected components of G located 
in the {0 < Im£, < Re£ < ir} are described in section 7.2.1, part 2. 

In Fig. 11, the curve Co intersects two connected components of G, the segment [(2^-1 ) C2«-2] of the 
imaginary axis and the segment [C2n>C2n+i] of the real line. So, Lemma 9.5 implies that 

m(c , 0, k 9 ) = m(c , Co, « 5 ) " 




Figure 11: The curve cq 



(9.19) 



7T 



Tn-lj 



ReAt 9 (C2n)) 



as ReK stays constant on any connected component of G. As n g is defined by the formulae (9.17) 
and (7.10), one has 

(9.20) «<?(C2ra) = ~«(C2n) = ~(«p(C2n) ~ TTTl) = - (7171- 7m) = 0. 

Along the interval [— C2n 5 C2n] ; one has re s (C) = ~ K (C) 6 ^! hence, 



« S (C2n-l) = -«(C2n-l) = -(«p(C2n-l) - = -(7T(n - 1) - 7Tn) = 7T. 
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Substituting this and (9.20) into (9.19), we finally get 

m(g, h,di) = m(c , 0, K g ) = 1. 
The result. Proposition 9.2, formulae (9.2) and (8.3) imply that, for £ G Si, 

(9.21) b = 6ie 2 ? i (H-o(l)), bi = exp (-^ J ndC, - + ^ fi_ + iA arg where /3 = /3 ,i- 
Completing the analysis. The coefficient &o being e-periodic, we write its Fourier series 

(9.22) 6 (C)= V 6 ,e 27r/C/e where b ,i = - b (C)e- 27Tl<:/e d( for Z G Z, 



As bo is analytic in the strip {|Im£| < Yq}, £ can be taken arbitrarily in the strip {|Im£| < Yq}. 
The asymptotics and the estimates for bo in Proposition 9.1 are obtained by analyzing its Fourier 
coefficients. To estimate the Fourier coefficients with non-positive index, one uses (9.18) and (9.22) 
with (" G Sq. To study the Fourier coefficients with positive index, one uses (9.21) and (9.22) with 
C G Si. We omit the elementary details and note only that bo in (9.18) is the leading term of the 
asymptotics of &o,Oi an d that b\ in (9.21) is the leading term for &o,i- 

9.6. The asymptotics of the coefficient a,Q. By (9.2), it suffices to compute the Wronskian 
w (fo(~7 C)> /tt(") C))- The computations of the coefficient ao follow the same scheme as the ones of 
bo- So, we only outline them. Now, 

(9.23) h = f:, g = f ; D h = V^ D g = V ; 

(9.24) Ch = vr, ( g = 0; 

(9.25) K h (C) = -7s(C) for ( ~ vr, and K g (() = -k(() for C ~ 0. 
Recall that the complex momentum is real on [C2n+i> 27r — C2n+i]- This imply that 
(9-26) K h (0 = -«(C) for C ~ vr. 

9.6.1. The asymptotics in the strip So. In this case, the meeting domain do is the subdomain of the 
strip So located between the lines the lines 71 and 72 symmetric to 72 with respect to the real line 
(see section 9.5.1). These two lines do not intersect. 

The arc 7(5, h, do) is the curve ao,o shown in Fig. 10. One has m(g, h, do) = 0. 
The asymptotics of ao for £ G So is described by 

(9.27) ao = a o(l + 00 = exp ^— -J nd( + J 0_ + iA arg g| a ^ where a = ao,o- 

9.6.2. The asymptotics in the strip S\. Now, the meeting domain d\ is the subdomain of the strip Si 
located between the line 73 (see section 9.5.2) and the line 72. 

The arc j(g,h,di) is the curve qo,i shown in Fig. 10. One has 

j(g,h,di) = co + j(g,h,d ), 
where cq is the closed curve shown in Fig. 11. The computation done for bo in Si yields 

m(g,h,di) = m(co,0,K g ) = 1. 
In result, for £ G Si, we get the asymptotic formula 

2Tvic (if 2iT 2 i f \ 

(9.28) ao = die = (l+o(l)), di = exp I — J Kd( h J 0_ + iA arg q\ a J where a = ao,i- 

The asymptotics (9.27) and (9.28) imply the formulae and the estimates for ao in Proposition 9.1. 

10. Phase integrals, tunneling coefficients and the iso-energy surface 

In this section, we first check the statements found in section 1.3.3. We also prove Lemma 9.3 
giving a geometric interpretation of the vertical tunneling coefficients. 

Then, we analyze the geometry of the iso-energy curves T and Tk (see (0.4) and (0.3)) and justify the 
interpretation of the phase integrals and tunneling coefficients in terms of these curves. 
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10.1. The complex momentum on the integration contours. The phase integrals and the tun- 
neling coefficients were defined as contour integrals of the complex momentum along the curves shown 
in Fig. 3 and 4. We have claimed that, on each of these curves, one can fix a continuous branch of the 
complex momentum, which we justify in 

Lemma 10.1. Let 7 be one of the curves 70, j n , 7/1,0 ; 7/1,71-; 7u,o an d %,tt- Any branch of the complex 
momentum, analytic in a neighborhood of a point of j, can be analytically continued to a single valued 
function on 7. 

Proof. The curve 7 goes exactly around two branch points of the complex momentum. They are 
of square root type (see section 7.1.1). So, it suffices to check that, at the branch points, the values 
of the complex momentum coincide. For the curve 7/1,70 this follows from the facts that £ (defined 
in (7.8)) bijectively maps the interval \(,2n-> C,2n+i] onto the n-th spectral gap of the periodic operator, 
and that the values of a branch of the Bloch quasi-momentum coincide at the ends of a gap. For 7^, 
this holds as £ maps the interval (C2n+i ; 2-7T — C2n+i) into the n-th spectral band so that both ends are 
mapped on E^n+i- For % j7r , it holds as £ maps the segment (C2n+2 ; 27r — C,2n+2) into the (n + l)-st 
spectral band so that both its ends are mapped on £^1+2- The analysis of the other curves is done in 
the same way. □ 

10.2. Independence of the tunneling coefficients and phase integrals on the branch of the 
complex momentum in their definitions. The independence follows from the observations: 

• only the signs of the integrals defining the phase integrals and the tunneling coefficients depend 
on the choice of the branches of the complex momentum being integrated; 

• the branches of the complex momentum being chosen, each of the phase integrals and each of 
the tunneling action is real and non-zero. 

Let us check the first observation. Let 7 be one of the curves 70, 7V, 7/^0 ; 7/i,7n lv,0 and % t7T . Let 
ac be a branch of the complex momentum continuous on 7. The formula (7.5) describes all the other 
branches continuous on 7. As 7 is closed, this shows that only the sign of the integral <f ndC, depends 
on the choice of the branch k. 

Recall that k p is analytic in the strip S p (see section 7.3). To prove the second observation, we fix a 
branch of the complex momentum on each of the integration contours. For 7^ and u , we fix this 
branch so that k = K p — irn on the parts of the contours in C+; for j vu , we choose K = k p — irn on 
the parts of the contours in C + fl {v < Re£}. We orient the contours %, 7/^ and % )7T clockwise, and 
we orient the contours 70, 7/1 an d %,o anticlockwise. Then, the second observation follows from 

Lemma 10.2. For E G J, for the above definitions of the integration contours and of the branches of 
the complex momentum defined on them, each of the functions & u , Sh,u an d S V)U is positive. 

Proof. Begin with & n . As C2n.+i is a square root branch point of k, and, as K(^2n+i) = 0, we get 

f27T — ^271+1 



r-ZTT — (,2n+l 

$AE)= / K(( + i0)d(, 



'C2n + 1 

where one integrates along R. As £(C) i s even, one proves that 

(10.1) $ n (E) = 2 K(( + i0)d( = 2 (K p (()-7rn)d( 

J C2n+1 •> C2n+1 



Inside the integration interval, one has ImK p = 0, and im < ReK p < 7r(n + 1). This implies the 
positivity of 

Arguing as above, for Sh,ir, we get 

/•C2n+1 

(10.2) S h>7r {E) = -i (ASp(C) - vrn) d(, 

where one integrates along R. Inside the integration interval, one has ReK p = irn and Iuik p > so 

that Sh,n > 0. 

For S v -n- , one obtains 



(10.3) S V)7T (E) = -2i / (Kp(C) - 7r(n + 1)) 

J(2n + 2 
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where one integrates along 7r+i]R. Inside the integration interval, irn < Re k p < 7r(n+l) and Im k p = 
which implies S V)% > 0. 

Arguing similarly, one proves the positivity of $o 5 S V) o and Sh,o- We omit further details. □ 

10.3. Proof of the inequalities (1.4). One has 

& n (E) = 2 («p(C) - vrn) d( and $ (E) = -2 / (k p (C) - vrn) df . 

•'C2n + l ■'0 

The first equality was established when proving Lemma 10.2. The second is proved similarly. In view 
of (7.9), we get 

/•7T rC,2n 

& n {E) = 2 k' p (E -ocosC)dC and ^o( E ) = -2 k' p (E - a cos Q d( , 

where fc p is the main branch of the Bloch quasi-momentum described in section 6.1.2. As, inside any 
spectral band of the periodic operator Hq, the derivative k' p is positive, this proves (1.4). 

10.4. Proof of (1.9). We can choose the oriented contours 7^0 and 7/^ so that one be the symmetric 
of the other with respect to the origin. As S (Q is even, for Q G 7/i j7r , one has «(— C) — K (0- These 
two remarks imply relations (1.9). 

10.5. Proof of Lemma 9.3. We shall prove the statement of Lemma 9.3 for v = tt. For v = the 
argument is similar. As S V<K (E*) G M, (10.3) implies that 

pit 

(10.4) S V , V (E*) = ReS v , n {E*) = 21m / (k p (() - vr(n + l))dC- 

Let us deform the integration contour in the right hand side so that it go successively 

• from C2n+2 along the Stokes line "b" to (b a , the point of intersection of the Stokes lines "b" 
and "a" (see Fig. 8), 

• from (b a along the Stokes line "a" to C2n+i, 

• from C2n+i to 7r along the interval [C2n+ii ^ which also is a Stokes line. 

As /%)(C2n+i) = ^ n and Kp{(,2n+2) = tt(w + 1), the definitions of the Stokes lines then imply that 

S V ,«(E*) = 21m / (ASp(C) - vr(n + 1)) d( + 21m / (k p (C) - ?r(n + 1)) d( 

J(2n+2, along "b" J(ba, along "a" 

+ 21m / (ASp(C) - 7r(n + 1)) = + 2^Im Cia + = 27rlm C fea - 

J&n+i, along R 

As the set of the Stokes lines is symmetric with respect to both the real line and the line ir + iK, the 
definition of Y v ^ implies that Im^ a = Y V ^(E*). This and the result of the last computation imply 
that 5^ )7r (£'*) = 2nY Vt n(E*). The proof of Lemma 9.3 is complete. 

10.6. The iso-energy curve. The iso-energy curve T is defined by (0.4). A point (£ k) € C 2 belongs 
to r if and only if k is one of the values of the complex momentum at the point £. 

We now discuss the iso-energy curve under the assumptions (H), (O) and (TIBM). 

10.6.1. The real branches. Consider the real iso-energy curve Tk defined by (0.3). Its connected 
components are the real branches of the iso-energy curve. One has 

Lemma 10.3. The real iso-energy curve is 2% -"periodic in both the k- and £- directions; it is symmetric 
with respect to each of the lines {k = nn} and {( = irm} for m, n £ Z. 

Any periodicity cell contains exactly two real branches of T. Each of them is homeomorphic to a 
circle. 

There exists 70 and j n , two disjoint connected components of such that the convex hull of 70 

contains the point (0,7rn), and the convex hull of 7^ contains the point (jr,irn). 

The curves 70 and 7^ are disjoint and are inside the strip {it{n — 1) < k < 7r(n + 1)}. 

Any other real branch of T can be obtained either from 70 or j n by 2it -translations in k- or/and in 

C- directions. 
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Proof. The analysis of T]r is quite standard. A detailed example can be found in [10]. So, we only 
outline the proof of Lemma 10.3. The periodicity and the symmetries of in £ follows from the 
symmetry and periodicity of the cosine and from formula (7.5). 

Describe two real branches of T. Recall that one has K p ([£2n+i> tt]) C [7m,7r(n + 1)[, /%>([0, C271]) C 
]7r(n— 1), 7m] and k p ([C2td C2n+i]) C 7m+zM + . On the first two intervals, K p is monotonously increasing; 
on the last interval, the imaginary part of k p has only one maximum; this maximum is non degenerate. 
The graphs of k p on each of the intervals [0, £20] and [£271+1, tt] belong to Tr. The real branch 70 is 
obtained from the graph on [0, ^2n] by the reflections with respect to the lines {k = nn} and {£ = 0}. 
The real branch 7^ is obtained from the graph on [C271+1 , 7r] by the reflections with respect to the lines 
{k = 7m} and {£ = tt}. 

We omit further elementary details of the proof. □ 

10.6.2. Complex loops. We prove 

Lemma 10.4. The closed curve 70 (resp. jh,o> 1h,-k, %,o and (see figures 3 and 4) is the 

projection on the Q-plane of a loop 70 (resp. j n , jh,o> lh,-w, 7^,o an d Jv,n) that is located on T. These 
loops satisfy: 

• the loop 7^ i7r connects the real branches "f n and 70; 

• the loop 'jhfi connects the real branches 70 and j n — (2-7r,0); 

• the loop 7„ i7r connects the real branches j n and 7^ + (0, 2tt); 

• the loop j V) q connects the real branches 70 and 70 + (0, 2n). 

In Fig. 2, we sketched the loops described in Lemma 10.4. 

Proof of Lemma 10. 4- By Lemma 10.1, the complex momentum can be analytically continued along 
each of the above closed curves on C. This implies that each of them is the projection to C of a loop 
on r. Fix v G {0, tt}. For d G {h, v }, the loops discussed in the lemma satisfy: 

(10.5) 7 „ = {(C, Kp(C)); C G %}, and ld , u = {(£ k p (Q)] ( G % u }. 

Here, for j VjV , k p denotes the branch of the complex momentum that coincides with k p on the parts of 
the contours in C + n {u < Re£}; for j u and jh,v, h is the branch that coincides with k p on the parts 
of the contours in C+. Therefore, we note that the curve 77^ intersects 70 and At the intersection 
point of jh,w arid 7^ (resp. 70), the branches of k p fixed on these curves coincide. This implies that 

connects the real branches 7^ and 70. 
The analysis of the other loops is done in the same way; we omit further details. □ 

10.6.3. Interpretation of the phase integrals and the tunneling coefficients in terms of the iso-energy 
curve. Let E be real. Pick v G {0,7r} and d G {v, h}. Formula (10.5) shows that, up to the sign, <fr u 
and Sd,v{E) coincide with ^ § nd( and — | § nd(. So, choosing the orientations of j v and in 

a suitable way, we get <E>j, = 5 § ndQ and Sd jV (E) = — | $ ndC,. 

11. Properties of the Fourier coefficients 

We now prove the estimates and the asymptotics of the Fourier coefficients found in Theorem 5.1 
which will complete the proof of this result. 

11.1. Computing the semi-classical factors. Proposition 9.1 shows that the leading terms of the 
first Fourier coefficients of a„ and b v contain factors of the form e~^~> Kd< ^. They are computed in 
Lemma 11.1. For E G J, one has 

(H.l) expf-1/ K dC )=e^t- h l 



£ 




■ <t> n -$,r-47r 2 (if I . ^Q+i-ir-iTT 2 



(11.2) expf--y Kd(\=e i 2, t v>0 t h ^, expf--^ Kd()=e i ^ ^,0^, 



(11.3) expl- [ 



(11.4) expfj/ nd( 
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Here, the integration contours are the curves shown in Fig. 10; in each of the integrals, k is the 
branch of the complex momentum obtained from the one introduced in the beginning of the section 7.3 
by analytic continuation along the integration contour from its beginning to its end. 

Proof. All the formulae (11.1) - (H-4) are proved similarly. Check the first formula in (11.1). 
Therefore, we deform the curve ao,o so that it go along the real line going around the branch points 
<^2n and (,2n+i along infinitesimally small circles. We get 

k d( = h + 1 2 + h 

"0,0 

where 

(11.5) h = - / fs(C + t'0)dC, h = ~ K(( + i0)d( and I 3 = - k((-i0)d(. 

Here, in I\ and I2, we integrate the branch of the complex momentum k introduced in the beginning of 
the section 7.3, and, in I3, k is the branch obtained from k by analytic continuation from the interval 
(C2n+i>7r) + iO to the interval (C2n+i> vr) — iO around the branch point C2n+i in the anti-clockwise 
direction. 

Consider I3. As (,2n+i is a square root branch point of k and as re(C2n+i) = 0, we have k(C — iO) = 
-k(( + iO) for C G (CjtH-i. C R. So, h = j k(C + iO) dC = / («p(C) - tm) dC- Comparing 

^ C2n + 1 J C2n+1 

this with the right hand side of (10.1), we get I3 = Similarly, one proves that I\ = ^3>o- I n view 

of (10.2), one has I2 = —iS^^. Combining the obtained expressions for ii, I2 and ^3, we get 

exp ^- % - J k d(j = exp (^(h +I 2 + h)^ = exp ^ ($„ + ^n) + ^ ^.Tr^ • 

This and the definition of th, n implies the first formula from (11.1). The second formula is proved 
similarly. 

Describe the computation of the integrals in (11.2). Let J ndC, be one of them. First, using a symmetry 
argument, we rewrite the integral in terms of the branch k p . As k is real analytic in a neighborhood 
of 0, one notes that KdQ = f-nd(, where 7 is the oriented contour symmetric to 7 with respect to 
the real line. One expresses the integral f^ndC, in terms of the tunneling actions and phase integrals 
using arguments similar the ones presented above, and, then one computes ndC, using the fact that 
the phase integrals and the actions are real for real E. We omit further details. 

Describe the computation of the integrals in (11.3) and (11.4). Let f ftd£ be one of them. Again, 
using a symmetry argument, we rewrite the integral in terms of the branch K p . As the function 
£ — ► K>(iC) is real analytic in a neighborhood of 0, one notes that J ndQ = —J_^ ndQ, where —7 is the 
oriented contour symmetric to 7 with respect to the imaginary axis. Then, one computes the integral 
J_-Kd<^ as the integrals in (11.1) and (11.2). We omit further details. This completes the proof of 
Lemma 11.1. □ 

11.2. Proof of (5.10) — (5.12). Being valid for E G J, formulae (11.1) - (H-4) remain valid in some 
neighborhood of J independent of e (as equalities between analytic functions). The formulae (5.10) 
and (5.12) follow from the asymptotics (9.5) and (9.6), and from formulae (11.1) - (11.4). To illustrate 
this, let us prove the formulae for oo,o- Let Vq be the neighborhood of E* from Proposition 9.1. 
Using (9.5) and (11.1), for E G Vq, we get 

(11.6) 

fl o,o = exp + $0) + j fi- + *AArg(/ Ko + o(l) \ = t^ Q exp (^(^ + $0) + 0(l)j , 

where we have used (1.10) and the fact that 0_ and q are independent of e. As E 1— > tf l ^ 7T (E), 
E 1— > $>q(E) and E 1— > $> n (E) are real analytic, (11.6) implies the representations concerning aoo 
from (5.10) and (5.11). 

11.3. Proof of (2.23). Pick v £ {0,7r}. Let V* be the neighborhood of E* from Theorem 5.1. By 
means of (5.44) (5.11) and (5.12), for E G V*, we get z v = 0(l/e). The Cauchy estimates then imply 
that z' v = 0(1/ e) in any fixed compact of V*. So, at expense of reducing somewhat V*, we have 
proved (2.23). 
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(a) 



(b) 



Figure 12: The curves g and g. 



12. Combinations of Fourier coefficients 



Here, we study the asymptotics of the quantities 0, Th, T Vj q, T V)T , <&q, ^>tt and zq, z w . We always use 
the branches k, ip± and 0± described in the beginning of section 9.2. Also, we systematically use the 
notations and constructions from section 6. 

Let So be a point in S. Assume that it is not a branch point, and that tt(£o) 6 K- Consider U, a 
neighborhood of £$ where 7r _1 is analytic. On U, we define the mapping * : £ \— ► 7r _1 (7r(£)). For 7, 
an oriented curve in S containing no branch points and beginning at £q, we continue the map * along 
7 and, thus, define the oriented curve 7*. 

12.1. The constant and the coefficients T^, T Vi o, T vir . They are defined in (5.41) and (5.42). 
The asymptotics (2.21) and (2.20) are obtained in the same way; so, we justify only the asymptotic 
for 9. 

The proof that, for sufficiently small e, in the case of Theorem 2.2, one has (2.21) with the constant 
9 n defined in (6.4), consists of three steps. 



12.1.1. Asymptotics of | ^ 
in Fig. 12, part a, and that, 



Let g be a curve on S that goes around the branch points as shown 
:*or 7r(£) > 0, is on the sheet of S where k p (ir(£)) is the Bloch quasi- 



momentum of ip(x,£). We check that 

a o,o 



(12.1) 



«7T,0 



exp / n(£) + / Q(£) + o(l 



g* 



The representations (9.5) and the first formulae from (11.1) and (11.3) imply that 



(12.2) 



«o,o 

«7T,0 




MO 



MO + 0(1) 



M0 + °(i) 



<*7T,0 



as t h ^ = th,o, see (1.9). 

Recall that the curves (cejy,o)i/e{o,7r} are shown in Fig. 10. We can and do assume that —0^,0 is the 
symmetric to ao,o with respect to the origin. 

As there are only two different branches of £ 1— ► 0(0) an d as the branch points of O coincide with 
those of k, the analytic continuation of + along c%o> near 0, the end of a n fi, coincides with 0_. 
Therefore, (12.2) can be rewritten in the form 



(12.3) 



1 a o.o I 




1 a-n-,0 1 





cxp 



0_ + f 0_ + o(l) 



Now, we make the change of variables £ 1— > 5(0- It maps each of the curves ao,o and — a^o on g, and 
we get 

ex P (/«o,o n -(0 + /-«.,„ n-(0) =exp (2j g n(£)) , 

where we have used that, for £ near 0, the branches £ — ► O±(0 correspond to the Bloch solutions 
^ i — ► ip±(x, £ (0) with the quasi-momenta C l—> =t^p( , ?(C))- I n section 6.2, we have formulated general 
properties of O. The fifth property implies that 



(12.4) 



0(£) 



0(f). 
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So, 

I ex P (Xo n -(0 + /-^ MC)) I = exp (j g m + n(£)) . 

This and (12.3) imply (12.1). 

12.1.2. Computation ofdetT n . Here, we prove that 

(12.5) det T n = — exp (- j VL{£) + Q(S)j . 

Relations (9.3), (8.3), and (8.4) imply that 

/ 19R \ r] . T _ »(/ T j;)l (+ 2, _ k' p (E+a)w(iP + ( - ,E+a),Tp_(- ,E+a)) 

[ 0) aect*- - k' p (E-a)w^+(-,E-a)^-(-,E-a))- 

Furthermore, it follows directly from the definition of O that Sl(£)+Sl(£) = — (flog Jq i/j(x, £)tp(x, £) dx. 

Note that ip(x,£)ijj(x,£) remains the same when we interchange £ and £. Therefore, it depends only 
on E = 7r(£) and is single valued on the complex plane. So, we get 

(12.7) exp (jn(£) + n(£)) = ^ + 



J9 ) J Q ip + (x,e)ip-(x,e)dx\ e=E+a ' 

On any simply connected domain of C containing no branch points of tp, one has (see, for example, [11]) 

l 

ifj + (x,E)ii_(x,E) dx = -ik'(E)w(^ + (-,E),^(;E)), 



where ip± are two different branches of ip and k is the Bloch quasi-momentum of "0+. This for- 
mula, (12.7) and (12.6) imply (12.5). 

12.1.3. Completing the proof of (2.21). Let g n C S be the curve shown in Fig. 12, part b; its part 
marked by "*" is on the part of S where k p {Ti{£)) is the Bloch quasi-momentum of ip(x,£). Rela- 
tions (12.1), (12.5) and (5.41) imply that 9 = exp (§- n tl(£) + o(l)) . 

Now, let us compare <f~ £l(£) with <f> g Q(£) where g n is the curve in (6.4). Note that, on S, modulo 
contractible curves, one has g n = g n . When deforming on S the curve g n to g n , one may intersect 
poles of f2. The poles and the residues of f2 are described in section 6.2. This description implies that 

the above two integrals coincide modulo 2iri. So, we have 6 = exp (^<f g Q(£) + o(l) \ . This completes 
the proof of (2.21). 

12.2. The phases {&v}v=o,n and {z u } u= q j7T . These are defined in (5.40) and (5.44). The asymptotics 
of all the phases (see (2.19) and (2.22)) are obtained in the same way; we justify only the asymptotic 
for 

So, we prove here that, for sufficiently small e, in the case of Theorem 2.2, $„- admits the asymp- 
totics (2.19). 

The asymptotics (9.5) and (9.6) and formulae (11.1) and (11.3) imply that 

(12.8) l$ w = lif> % + ±(s-S) + ^8 + o(l), 
where 

(12.9) s=f n++f n_+r, n+-L n_, 

( 12 - 10 ) s = Aargg|^ o + Aargg| aoo + AarggL^ - Aargg^, 

where (a Uj o, Pv,o) e e{o,ir} are sketched in Fig. 10. We can and, below, we assume that, as the oriented 
curve «o,o (resp. /?o,o) is symmetric to the oriented curve — a n< o (resp. — /^o) with respect to zero. 
First, show that S — S = 0. Arguing as when deducing (12.3) from (12.2), we get 

Now, we make the change of variables £ i— > £(()■ As £(ao,o) = £(—&n,o), and £(—P n ,o) = £(A),o)> we 
get 
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As when proving (12.4), we see that the terms in (12.11) differ only by complex conjugation. So, 5 is 
real and 5 — 5 = 0. 

Finally, we show that that s = 0. This will complete the proof of the asymptotics of ^ 7T . 
When computing the increments of the argument of (/(£) = y/k' (£((,), we choose the (continuous) 
branch of this function which is positive on the interval (— C2niC2n)- Then, in a neighborhood of 
zero, q*(C) = q(C) an d <?(— C) = q(0- Therefore, and due to our "symmetric" choice of the curves 
(av,oW{o,7r} and (^,o) we { ^}, we get 

Aargg|^ o = - Aargq\ aQ0 and A arg^ Q = - A arg q\^ = A arg q\ po fi . 
This and the definition of s implies that s = 0. □ 
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